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Abstract 

We consider the stationary Stokes problem in a three-dimensional fluid domain T with 
non-homogeneous Dirichlet boundary conditions. We assume that this fluid domain is the 
complement of a bounded obstacle B in a bounded or an exterior smooth container St. We 
compute sharp asymptotics of the solution to the Stokes problem when the distance between 
the obstacle and the container boundary is small. 
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In this paper, we consider the 3D-Stokes problem 

Au-Vp = 0, (1) 

Vu = 0, (2) 

in a fluid domain T. Without restricting the generality, we set the viscosity of the fluid to 1. We 
assume that T = SI \ B is the complement of a smooth simply-connected bounded domain B inside 
a container SI. The container SI is either a relatively compact simply-connected smooth open set or 
the exterior of a simply-connected smooth compact set B*. In both cases, it has a smooth compact 
connected boundary. We complete then (ITjl -(f2 jl with boundary conditions: 

u = u*, on dSt , (3) 

u = 0, on dB , (4) 

where u* £ C°°{dS I) does not prescribe any flux through the container boundary: 

[ u* • nder = 0 . (5) 

Jan 

If SI is an exterior domain, we add a vanishing condition at infinity: 

lint up, y,z)= 0. (6) 

\{x,y,z) |- kx > 

With the above assumptions on the boundary data u* it is classical that system © ©(+©) 
admits a unique classical solution (u,p) (the pressure p being unique up to a constant), see [8] for 
instance. Our aim in this paper is to give a sharp description of the solution to the Stokes problem 
© ©(+©) when the distance between B and dSl is small. 

Computing such asymptotics is an important issue related to the modeling of solid-body motion 
inside a viscous fluid. The typical configuration we have in mind is that B (resp. B and B*) is a 
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(are) moving solid body (bodies) inside the container Q (in the whole space) which is filled by a 
viscous incompressible constant-density fluid. A typical issue is then to determine whether the fluid 
viscosity prevents the moving body B from touching other solid boundaries ( i.e. the boundary 
of the container or the boundary of the other solid body) and more generally to measure the 
influence of the viscosity on the close-contact dynamics of B. To this end, one remarks that, when 
a solid body is about to collide another solid boundary with moderate relative velocity, the fluid 
Reynolds number tends to 0 in the gap so that a stationary Stokes system is sufficient to predict 
the force and torque exherted by the fluid on the moving body. With this particular application 
in mind, several authors consider the free-fall of a sphere above a ramp H El El Hi DU in a Stokes 
fluid. Explicit values for the solution to the Stokes problem £[]) (J4]) (-T dHI) ) and the associated 
force and torque are provided. A formal lubrication approximation is also proposed in mm 
which generalizes these formulas to arbitrary configurations. In the limit regime where there is 
contact, solutions to the Stokes problem are also computed in m under further assumptions on 
the boundary data u* (broadly, the boundary data u* has to vanish sufficiently where there is 
contact). Related computations for a perfect fluid are provided in [1,, 15, IB] . 

In this paper, we fill the gap between the explicit formulas of ia n 0 na ITB] , the formal 
asymptotics of mum and the analysis in m- We justify rigorously lubrication approximation 
in the spirit of M- Compared to these latter references, we are interested inhere in fluid films 
that do not vanish uniformly in their widths. This fact leads to severe new difficulties. First, 
the lubrication scaling acts on coordinates in both tangential and orthogonal directions to the 
boundaries (see © (ED)- Second, the asymptotic pressure and velocity-fields yielding from the 
formal lubrication approximation are defined on an asymptotic fluid domain which is not simply 
related to the fluid-domain for a given positive body/boundary distance h. Consequently, in order 
to compare the values of the solution to the Stokes problem with its asymptotic value, we introduce 
an intermediate velocity-field which embeds the lubrication approximation in the effective fluid- 
domain (see Section |3.1l) . The construction of this intermediate velocity-field is an important 
step of our analysis. Indeed, the intermediate velocity-field is a key-ingredient in order to extend 
the computations on the close-contact dynamics of bodies in a Stokes fluid to more complicated 
models: Navier-Stokes/Newton models or Navier-Stokes/elasticity models. This fact has already 
been shown for the Navier-Stokes/Newton model in simple configurations (see [SJCE2] 1131 EH!)• We 
believe our approach extends to other incompresible models (for instance to the case of potential 
flows as in jlB]). However, the incompressibility condition is crucial to our computation and the 
extension to the compressible case is completly open. 

To fix ideas, we make more specific the geometry of the gap between B and <9fh We introduce 
a set of cartesian coordinates ( x , y, z) G R 3 and (r, 6, z) G (0, oo) x (— 7 r, 7 r) x R the corresponding 
cylindrical coordinates. These coordinates are associated with two orthonormal basis of R 3 denoted 
by (e x ,e y ,e z ) and (e r ,eg,e 2 ) respectively. We consider that, in a neighborhood U of the origin, 
Q and B satisfy: 

{x, y, z) G (R 3 \ Cl) n U ^ {(x, y) G 5(0, L) and z < 7 b(x, y)} , (7) 

(x,y,z)GB DW {(x,y) <E 5(0,5) and 2 >/i + 7 t (x,?/)} . ( 8 ) 

Here (5, L) are given positive parameters. The first one measures the distance between dfl and 
dB, the second one is a characteristic length on which parametrizing the boundaries of Q and B 
by (x, y)-variables is relevant. The two functions 7 1 and 75 are smooth on 5(0, L) C M 2 and we 
assume throughout the paper that they satisfy the following assumptions: 

(Al) V 7t (0) = V 7b (0) = 0, 

(A2) there exists i?i > 0 such that (in the sense of symmetric matrices) 

V 2 7 t(a:, y) ~ V 2 7 b(x, y) > 2-I 3 , V (x, y) G 5(0, L). 
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We complement these local assumptions in U with a global one concerning B and <912 outside U : 

3J>0s.t. dist(dfl\lA,B) > S. (A3) 

The parametrizations (0 and © together with assumptions (Al)-(A2)-(A3) are paradigmatic 
of what we call a ’’single non-degenerate contact”. Indeed, considering that B is a translating 
particle inside the container 11 amounts to let the parameter h depend on time. Assumption (A3) 
implies then that a contact between B and <9f1 may only hold in U. Assumption (A2) yields that 
the contact in IA is unique and non-degenerate i.e. when h = 0 the ’’vertical” distance 7 1 — 7 b 
vanishes in 0 only and with minimal vanishing order. We emphasize that, in the smooth case we 
consider here, the uniform boundedness we require in (A2) reduces to assuming that 

V 2 7 t(0) - V 2 7b (0) > i-Ia , (A’2) 

Ki 

up to restrict the size of L and change the value of R\. In this ’’single non-degenerate contact”- 
case, the assumptions (|7 ]l -([8 ]> together with (Al) do not restrict the generality: they only amount 
to choose the origin of coordinates in the only point of <912 realizing the distance between 1 and 
<9 B, and to choose the system of coordinates so that the common normal to 512 and dB, in the 
pair of points realizing the distance between <912 and dB , is e z . 

Our aim is to compute the asymptotics of the solution to the Stokes problem © ©(+©) 
for a given boundary condition u* e C'°°(911) in the geometry depicted above under the further 
assumption that h is small and other parameters are of order 1. To introduce our main result, we 
recall the main steps of the formal computations in for the case where B is a sphere of radius 
S and = M. 3 \B* with B* a sphere of radius R. First, given the shape of the aperture between 
both spheres, one looks for a solution (u,p) := ((u x , u yi u z ),p) that reads 


u x (x,y, z) 

= h a ~^u x (h~^x 1 h~^y,h~ 1 z), 

(9) 

Uy(x,y,z) 

= h a ~^u y (h~^x,h~^y,h~ 1 z), 

(10) 

u z (x, y, z) 

= h a u z (h ~2 x,h~2 y,h~ 3 z ), 

(11) 

p(x, y, z) 

= h a ~ 2 p(h~^x,h~^y,h~ 3 z), 

(12) 


in the aperture between the spheres. The parameter a is chosen depending on the values of 
u* := (u*,u*,u*). For instance, if u* = u*\_e z with u* ± £ R \ {0}, one chooses a = 0. We proceed 
with this particular case. We denote with tildas the new space variables 

x := (h~2x, h~^y, h~ 1 z). 

These new coordinates belong to the set Qh that ’’converges” (when h — > 0) to 


glub | (x,y,z) € R 3 s.t. 5 S 
Substituting ansatz © dll into m m yields: 


x 2 + y 2 x 2 + y 2 \ 2 
2 R ’ + 2 S )]' 


dx “f dy Uy T d z u z — 0, 

dzz'&x dxp = 0 , 

dz Z Uy d y p — 0 , 

dzp = 0 , 

completed with boundary conditions: 


Si 

8 

II 

Si 

II 

0 

on dg lub , 



0 

II 

N 

iS 

on dg l 2 ub ■- j 

|^ = 72 (x,y) =: 

i + * 2 2 + /.(-.a^ 2 } 

u z = U*J_ 

on dg[ ub ■- j 

= Ji(x,y) =: 
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The pressure is normalized by assuming that it vanishes at infinity. Introducing 7 := 72 — 71 , the 
unique solution to this problem reads u(x) = u“ b (x), pfk) = p lub (x,y) with: 


,lub 


{x,y,z) = ^d £ p lub (x,y)(z-j!{x,y))(z-j 2 (x,y)) 


u l y b (x, y, z) = -dyp lub (x,y){z-j 1 (x : y))(z^ ; y 2 (x,y )) 


u l ? b (x,y,z) = -div 2 iS 


[ ((S - 72(x, y))(s - 71 (x, y))) S7s,yP lub (x, y)ds 

J Z 


and where p lub is the unique solution to 

- ldiv(7 3 Vp z “ b ) 
lim p lub (x,y) 

\(x,y)\->oo 


u*l , onK 2 , 

0 . 


(13) 

(14) 


Herein, we justify these formal computations rigorously. A (nonetheless formal) statement of 
our main result reads: 


Theorem 1 . Assume (A1)-(A3) are in force and u* £ C°°(dQ) satisfies Let v* and be 
given by: 

v*/= u* x (0)e x + u* y (0)e y , v* ± (x, y) = u*(0) + xd x u* z (0) + yd y u* z (0 ). (15) 

When h « 1 while ( 7 t, 7 b, i?i, u*) remain of order 1, the main contribution to the velocity-field 
of the solution (u,p) to the Stokes problem (jTJ) (131) (+ (El) ) is given by v = (v//,v±) with 


v//{x,y,z) 


v±{x,y,z) 


t;(z - {h + 'Yt(x,y)))(z-7b(x,y))Vx, y q{x,y) + ^ + 11 , x ’ ^ — —v* 
2 j(x,y) 


v x , y 


h+ 7t (x,y) 


+ 


h+j t (x,y) 


div 


((s - 7 b (x, y))(s ~(h + 7 t(x, y)))) V x>y q(x, y)ds 
(h + 7 t(x, y) - s) 


x,y 


7 (x,y) 


d s 


(16) 

(17) 


in the aperture domain Gl /2 := {(x,y,z) £ R 3 s.t. ( x,y ) £ 5(0,L/ 2 ) and z £ ( 7 b(x,y),h + 
"ft(x,y))}, and where q is the unique solution to 


- — div (7 3 Vg) 


Re - ^div Xi j / (( 7 t + 7 b)v*), on B{ 0, L ), 

0 on (0, L ), 


(18) 

(19) 


with 7 := h + 7 t — 7 b- 

We introduce here the notations // and _L for the components of a vector v that are respectively 
parallel and normal to the tangent space to 9f l in the origin. Corresponding decomposition of V 
are denoted by V x<y and d z . We keep this convention in what follows. Our geometric assumptions 
imply that e x and e y are tangent to dQ in the origin. Hence, the derivatives 3 x u*(0) and 9 y u*(0) 
are well-defined. Even though the result concerns an asymptotic behavior when h —► 0 the value 
of this parameter is fixed in all computations. For this reason, we do not let the parameter h 
appear in most of our notations (such as the fluid domain _F, the distance function 7 ...). 

A more quantitative statement of the above theorem is given in next section. In particular, 
we make precise for which norms the extracted contribution is the dominating term, the size of 
remainder terms and the dependencies w.r.t. 7 *, 7 b, R\, u*. We prove this result by 

• computing a priori estimates on the pressure q as constructed in this statement (in the 
regime h « 1 ) 
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• relating the difference between the exact solution u to the Stokes problem and (an extension 
of) v to the obtained estimates on the pressure q. 

As a corollary, we obtain that the leading order in the asymptotics of the Stokes solution is given by 
the first order expansion of the boundary data u* in the origin. This dependance occurs through 
a pressure solution to a simpler problem but still not explicit. We explain in next section that in 
some special case (when B is a sphere and Q = R 3 \ B* with B* a sphere) we can compute accurate 
informations on the pressure and get explicit expansion of the Stokes solution w.r.t. h and the 
boundary data u,* (0), u)j_(0), V^u^O). 

In the set of assumptions we introduced up to now, we enforced the boundary condition to 
vanish on dB. We can always reduce asymptotic computations with general Dirichlet boundary 
conditions to this case. We also emphasize that the computations, that we present here in the 
’’single non-denegerate contact” case, extend to general ’’non-degenerate contacts”. A general 
non-degenerate contact would correspond to the case where the moving solid B may collide <90 
in several points satisfying assumption (A2). There would then exist at most a finite number 
of ” contact points” that might be treated separately as ” non-degenerate isolated contacts”. The 
linearity of the Stokes problem ensures that the exact solution behaves as the sum of the asymptotic 
solutions that we compute in the vicinity of each contact point. Previous calculations due to V. 
Starovoitov show that these ” non-degenerate isolated contacts” are the only remaining ones to rule 
out in order to show that no contact between solid bodies occur in a Stokes or a Navier-Stokes 
fluid (see [2D]). 

The outline of the paper is as follows. In next section, we recall the classical theory for 
solving the Stokes problem Jl]) |4])(-(-([6])) and give a quantitative statement of our main results, 
see Theorem 0] and Theorem [5] In Section [2] we study the properties of the problem (ITSll - dTDl) . We 
compute estimates satisfied by the solution q with respect to the data (u*, v* ) and compute the 
divergence-rate of q when h 0. One particular feature of the estimates we obtain is that they 
are ’’uniform” w.r.t. the distance h. The last two sections are devoted to the proofs of Theorem 
[4] and Theorem [5] respectively. 


1 Quantitative statement of main results 

In this section, we recall function-spaces and existence results for problem (JT]) (J3I) (+dHl)) • We give 
then a quantitative statement of our main result. We conclude by exhibiting two criterions which 
measure whether a velocity-field is a good approximation of a solution to the Stokes problem or 
not. 


1.1 Solving the Stokes problem 

We first recall the way the system ([T]) H])(-(-(l6|)) is tackled in [8|. Let denote: 

• V := (u^,u G C^QR 3 ) s.t. V • u = 0} and V 0 := {u € C™(F) s.t. V ■ u = 0} , 

• V (resp. Vo) the completion of V (resp. Vo) endowed with the norm: 


l|u:V|| 



Ur J 

This makes (V, || • ; V||) (resp. (Vo, || • ; V||)) to be a Hilbert space endowed with the scalar product: 


((u,v))= jf Vu: Vv. 

In the terminology of j8], our space V (resp. Vo) is a closed subspace of Z? 1 ’ 2 (J r ) (resp. Dq' 2 (jF)), 
see 0 p.80]. In particular, applying [H Theorem II.6.1 (i), (II.6.22)] we have that V embeds in 
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L 6 (J r ) so that V C W^(T). This entails that v £ V vanishes at infinity in a weak sense, if 
required, and has a well-defined trace V| a:F £ H^(dT). In particular, for arbitrary u* £ 
we might define the affine-subspace of V: 

C[u*] := {u £ V such that U| ae = 0 and U| aBi = u*} . 

We recall that this set is not empty as soon as u* prescribes no flux through <9f l and that we have 
the identity V[0] = Vo (see [5] Theorem II.7.1 (i)]). 

Following 0 we introduce the definition of generalized solution to (jTJl (J4J) (-(- (J6]) ) : 

Definition 2. Let u * £ H? (dfl) we call generalized solution to (jT]) ([i|l(+(p|)) any u £ V[u*] such 
that 

J Vu : Vv = 0, Vv £ Vo. (20) 

and we have the classical theorem (see jSJ Theorem V.2.1 and Theorem V.1.1]): 

Theorem 3. Given u* £ H^{dLl), such that : 


/ u* • nclcr = 0 , (21) 

Jd n 

there exists a unique generalized solution u to |l|) j4|)(-l-([6])). Moreover, 

• there exists p £ L^ oc (iF) such that (JT|) holds in the sense of distributions, 

• if u* £ C°°(dLl) there holds (u,p) £ C' 00 (.F; R 4 ) . 

A consequence of this result is that generalized solutions with smooth data are classical so¬ 
lutions. As we only consider this particular case throughout the paper, we drop the adjective 
’’generalized” in what follows. Also, we abusively call solution to ([l|) (0])(+([6])) a velocity-field u, 
even though a solution to the Stokes problem is a pair velocity/pressure (u,p). 


1.2 Main results 

The main contribution of this paper is the following theorem: 

Theorem 4. Assume (Al)-(A3) and boundary condition u* £ C°°(dQ) satisfies (1211) . Let denote 
by u £ V the unique associated solution to © ©(+©)• If h < 1, there exists v £ V and a 
constant K := K(Ri,Cb,S,dLl) such that: 

• v is given by (fTBl) - ([T71) in Gl /2 

. If || 7t ; C 3 (I?(0, L))|| + || 7b ; C 3 (B( 0, L))|| < C b and ||u*; H 3 (B(0, L))\\ + ||u* ; H* (9n)|| < C b 
there holds: 

ll u “ v;V|| < K (l + K(0)||ln(/i)|^ . 

If the distance function 7 is moreover radial, we have the following improvement of the second 
assertion: 

||u- v;Vj| < K. 

We remark that u and v are bounded uniformly in h as along as this parameter ranges a 
compact subset of (0,1]. Consequently, this theorem brings relevant informations in the limit 
h —> 0. It shows in particular that even though the fluid is incompressible (so that the Stokes 
problem is non-local), when h goes to 0 the leading order of the velocity-field is completely fixed 
by the boundary conditions and the geometric properties of the boundaries in the aperture Gl- 
It will be clear from the example given below that the velocity-field v we exhibit is larger than 
the remainder in the sense of the F-norm. We mention that our first interest in this problem was 
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to compute the diverging term in the case of spheres. As the remainder is bounded in this radial 
case, we did not look for a more complicated expansion of the solution. Nevertheless, a corollary 
of our method is the construction of the linear problem on which could rely the computation of a 
full expansion of the solution in terms of h. However, we would only describe the solution in the 
gap Ql and the computation of 0(l)-terms would require another tool enabling to compute the 
expansion of the solution outside the gap also. 

To highlight the relevance of the above theorem, we detail the case where B is a sphere and 
n = R 3 \ B* with B* a sphere, as in the computations in the introduction. Keeping the convention 
that B* has radius R and B has radius S, we have that, close to the origin, the functions jt and 
7 b satisfy: 

7 t(x,y) = X 2 s + Q (( a;2 + ^ 2 )^)> 

7 b(x,y) = ~ X ^R +°(( a;2 +y 2 )^), 


and also that 7 satisfies: 

7 (x,y) 


x 2 + y 2 
2Ri 


(x 2 +y 2 ) 2 
8 iZ| 


+ 0 ((x 2 + y 2 ) 3 ) 


with (i?i, R 3 ) g (0, oo) 2 given by: 


1 _ 1 1 11 1 
R^~S + R’ _ 53 ' 


( 22 ) 


In this case of spheres, we obtain: 

Theorem 5. Given a boundary condition u* £ C°°(dfl) satisfying (1211) there exists a constant 
K[R,S] depending only on R, S so that if h < 1, the unique solution u £ V to the associated 
Stokes problem (IB ©(+©), satisfies: 


ll u ; v\\ 2 


6tt|u1(0)| 2 



161?! 8 R 3 3Rf\ 
— -US-- 


+ A^lu^O)! 2 + ^R.IR.V^uKO) + iA-A u ;(0)| 2 l |ta(A)| 
+ ^[i?, 5 ](||u*;?? 3 (i?( 0 ,L))|| 2 + \\u*;Hi : dn\\ 2 ). 


The norms that we compute in this theorem are related to the forces and torques exerted by 
the fluid flow on the spheres. Indeed, assume that the sphere B is moving with a rigid velocity and 
that we want to compute the forces and torques exerted by a Stokes flow on B when the distance to 
the other sphere B* is small. By symmetry, we can assume that B* is moving with a rigid velocity 
u* and compute the force and torque exerted on B* . Then, we recall that the set of rigid velocities 
is a 6-dimensional vector space whose elements are characterized by a translation and angular 
velocity (computed with respect to the origin for simplicity). Let denote by (U-t^P-),;),-!;^ 
(resp. (Uo,i,Po,i)t=i,2,3) the solutions to the Stokes problem on T with elementary translational 
(resp. rotational) boundary conditions 



e, ; , on dB* , ( U 0 ,i 

0 , on dB, ( reSp - U 0 - 


ejXx, on dB *, \ 
0 , on dB , J 


for i = 1,2,3. 


Due to the linearity of the Stokes problem ([]]) d 4 ]) (+ (J6J) ), the solution (u,p) associated with 
boundary condition u* is then a combination of the (U->,j,-P->.,i)i=i,2,3 and (U07, -Po,*)»=i,2,3- 
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Furthermore, straightforward integration by parts yield that the j th component of the force exerted 
by the flow on B* reads: 


Fj := f (2D(u) — pl 3 )i/da ■ ej = 2 f D(u) : 
JOB- Jt 


Vu : VU_ 


IT 


(here v stands for the normal to dB* pointing towards B*). If u* = ej. we obtain: 


Fj = 



|VU_ 


and the asymptotic behavior or Fj when h —> 0 yields by applying the above theorem. By standard 
algebraic formulas this identity is generalized to all components of the forces and torques associated 
with any rigid velocity u*. Then, one can see the asymptotic expansions we compute in Theorem 
[5] as a justification and an improvement of the asymptotic values for the matrix K. provided in jBJ 
Section 7] (see [7.6]). 


1.3 Two approximation criterions 

We conclude this section by providing two criterions which will enable us to measure the distance 
between an approximation v and the exact solution u of the Stokes problem (JT]) (|l |) (+(p |) ). 

First, we recall that one way to prove the existence part of Theorem [3] is to construct u My £ 
V[u*] and to remark that we have V[u*] = u bdy + Vo- Existence and uniqueness of a generalized 
solution then yields from an application of the Stampacchia theorem. This proof entails the 
expected result together with the following variational characterization: 

Proposition 6. Given u* £ C°°(dfl) satisfying (l2lT) . the solution u to © ©(+©) realizes: 

||u;P|| 2 = min|y |Vv| 2 , v £ V[u*] 

As a consequence, given u* £ and u £ V the generalized solution to © ©(+©), 

we have conversely that any v £ P[u*] satisfies: 

IIv ; Vl| > ||u ; PH . (23) 

This enables to compute many bounds from above for ||u; V|| by choosing approximations v. The 
more relevant the approximation, the sharper the bound. To control the distance between an 
approximation v and the exact solution u we have more precisely: 

Proposition 7. Let u* £ C°°( dSY) satisfy (l2lT) and denote by u £ P the solution to © (|4|)(+([6|l). 
Given (v, q) £ P[u*] x C°°(F) uie denote: 

C\y, q] := ||Av-V g ;T 0 *||, 

i.e. C[v, q] is the best constant C such that: 

|x>'(^)((Av- Vq) ,w)- d( jf ) | < Cj|w; P|| , Vw £ V 0 . (24) 

Then, there holds: 

ll v — u ; V\\ <C[v,g], 

Proof. The proof is standard. For completeness, we recall it briefly. There holds: 



-u;P|| 2 = J |Vv — Vu| 2 


IT 


Vv : V(v-u)- J Vu: V(v-u). 


( 25 ) 


Applying (EU1) to w := v — u £ Vo, we obtain: 


Vu : V(v — u) = 0 . 


it 


Introducing then a sequence w„ £ Vo such that limn^oo w„ = w in V, we have: 


|| v — u; V\\ 2 = lim [ 

n—> OO J jr 


Vv : Vw„ 


where, for all n £ N, there holds: 


J Vv : Vw„ = J (Vv - gl 3 ) : Vw n = -((Av - Vg), w„) . 


Again, as w n £ Vo, we apply definition (l2dl) of C[v, q}. This entails: 


Vv : Vw, 


it 


< C[v,q ] ||w„; V\\ , 


and, in the limit n —> oo : 


||v-u;V|| 2 <C[v,g]||v-u;V||. 


This ends the proof. □ 

We emphasize that in the statement of this criterion there is no geometrical constant in front of 
C[v, q] in (12511 . This is particularly important as we want to consider the influence of the geometry 
on the relevance of the approximation v. Actually, these geometrical dependencies are hidden in 
the computation of C[v, g]. 


2 Preliminary results on the lubrication problem 


In this section, we consider the two-dimensional divergence problem: 


div [ 7 3 Vot] = /, on B( 0 ,i), 

(26) 

vj (x,y) = 0, on 913(0, L). 

(27) 


We restrict to source terms / £ C°°(B(0,L)) which have the special form: 

/ = w*-idiv(( 7t + 7 b )v*), where (v*,w*) £ C°°(B(0, L); R 2 x R). (28) 


The weight 7 £ C°°(B(0,L)) is computed with respect to 7t and 7 ^: 

7 ( 2 ;, y) = h + j t (x, y ) - j b (x, y ), V (x, y) £ B(0, L), 

with a small but positive distance h. The assumptions (Al)-(A3) are in force. It is standard that, 
since 7 is smooth and bounded from below by a strictly positive constant, there exists a unique 
smooth solution vj to m-m- Our aim is to compute estimates on quantities of the form: 


\S7 k w\ 2 . 


n £ N,fc £ N. 


JB(0,L) 

with explicit dependencies in the data /, the distance /i, 7 1 and "f b - In these estimates will appear 
constants depending directly on 7 , 7 t, 7 h- We state the definition of these constants as a lemma: 
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Lemma 8 . There exist constants ( C cvx , C e u) £ (0, oo) depending on 7 1 and 7 & only and (Cj^UeN £ 
(0,oo) N such that there holds: 


|| 7t ; C k {B{ 0, L))|| + ||7i; C k (B{ 0 , L))|| < C r k eg , 

C cvx < Aj(x, y), 


\/k£ N, 

V (x, y) £ 5(0, L), 


(29) 

(30) 


and, for h < 1 : 


h + C eU {x 2 +y 2 ) < 7 (x,y), V (x, y) £ 5(0, L) . (31) 

The proof of this lemma is an obvious consequence of (Al)-(A3) and is left to the reader. We 
remark that the constants C cvx and C e u are related to the constant R± appearing in assumption 
(A2). Assumption (Al) together with (Effl) imply there exists a constant K depending on C e u and 
C™ 9 such that: 


\^Tt{x,y)\ + \\7yb(x,y)\ < K\~j{x,y)\* , V (x, y) £ 5(0, L) , (32) 

Tt(x,y) + 7 b {x,y) < K^{x,y), V (x, y) £ 5(0, L) . (33) 

With these conventions, the main results of this section are the following propositions. We first 
write two weighted first-order estimates on the solution w of G51) - (l271) depending on the behavior 
of the data v* and w* (recall that / is given by (1281) 1 in the origin. 

Proposition 9. Assume that the data v*,w* satisfy: 


v*(0) = 0 , w*( 0) = 0, Vw*(0) = 0. 

Let n £ [0,oo). There exists constant K n depending on C cvx , C e u , C r 2 eg and n for which: 



|V^| 2 < K n { ||v*; 5 2 (5(0, L))|| 2 + |K; 5 3 (5(0, L))\\ 2 } 


(34) 


(35) 


Proposition 10. Let n £ [0,oo). There exists constant K n depending on C cvx ,C e u,C 2 eg and n 
for which: 

• if n = 0 we have: 


f 7 3+ "|vh 2 < ^{ll v *;^ 2 ( s (o^))ll 2 + lk*;^ 2 (^(o,A))|| 2 }, (36) 

Jb(o,l) n 

• if n £ (0,1) there holds: 

j l ) 73+ " |Vro|2 - Kn + llv*; H 2 (B( 0 ,L))W 2 + |K;5 2 (5(0, L))f|, (37) 

• if n = 1 we have: 

[ ^IS/w] 2 < K 1 {\w*(0)\ 2 \ln(h)\ + \\v*-,L 2 (B(0,L))\\ 2 + \\w*-,H 2 (B(0,L))\\ 2 } (38) 

JB(0,L) 

• if n > 1, there holds: 



l7| 3+ ”|Vti7| 2 

,L) 


<^„{||v*;L 2 (i?( 0 ,L ))|| 2 + ||u;*;77 1 (i3(0,L))|| 2 } 


(39) 
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We emphasize that the constants K n do not depend on h. In particular, the latter proposition 
implies that, when w*( 0 ) vanishes, all quantities: 



|7| 3+ ”|V?Z7| 2 , 

+) 


n £ N, 


remain bounded when h goes to 0. We complement the study with higher order estimates away 
and around the singularity point ( x,y ) = ( 0 , 0 ): 


Proposition 11. Given a integer k > 0 and 0 < e < L, there exists a positive constant K k es 
depending only on L, e, fc, C e u, C cvx , C k f 9 , with k! = max(fc + 1,2), for which m satisfies: 


\\uj; H k+2 (B(0, L) \ !?(0, e))|| < K k eg [||v*; H k+1 (B(0, L))|| + ||«/ 


■,H k '-\B(0,L))\\ 


(40) 


Proposition 12. Let k £ {0,..., 6 }, n £ (k,oo) and denote e*, := max(fc,3). There exists a 
constant Kff'f 9 depending on fc, n, C e u , C cvx , C k f 9 , with k' = max(fc + 1, 2), such that: 


[ | 7 | 3 +n| V k + 1 1 - 7 |2 < K™ 9 { f | 7 | 3 +(n-fc)| Vn7 | 2 + ||v*;lf fc+1 (B(0,L))|| 2 

Jb(0,L) ’ *- Jb(0,L ) 

+ \\w*-,H e '°(B(0,L))f}. (41) 

In the remainder of this section, we first give proofs for these propositions. We then conclude 
by showing that the first order estimates are optimal in the case of sphere. 


2.1 First order estimates : proofs of propositions l9l and 1101 

System (|26|) - (j27]l is associated with the weak formulation: 

^ / 7 3 Vn7 -V+= / ftp, V^C“(B(0,I)). (42) 

JB(0,L) JB(0,L) 

Hence, to construct solutions, one introduces the bilinear form: 

(<+V0 '-t ((+,^)) 7 : = +T / 7 3 V+-V^. 

Jb(o,l) 

A definition for weak solution to m-m then reads: 

Definition 13. We call weak solution to (tMI - dUD any w £ Hq(B(0, L)) such that: 

fv, V v? £ Hq(B(0, L)). 

Jb(o,l) 

Thanks to (P2!)l) - (l3Tl) . ((•, -)) 7 defines a coercive continuous bilinear form on Hq(B(0, L)). Ex¬ 
istence and uniqueness of a weak solution to (lMll - (E71l yields as a straightforward application of 
the Stampacchia theorem. Given our regularity assumptions on the distance function 7 , classical 
results on divergence problems apply so that this weak solution is smooth on B( 0, L) and satisfies 
(l26l) - (l27l) in a classical sense (see (TOj Chapter 8 ]). All the estimates in [TO] Chapter 8 ] depend a 
priori deeply on h so that an alternative approach is required. 

In the proofs below, we denote with K a constant which is important to our computations and 
shall put in brackets its relevant parameters (such as (C cvx , C e u), (C k e9 )k£f$)- These constants 
may differ from line to line. Most of these constants shall also depend on the parameter L but 
we include this tacitly. Notations C are used for generic constants that depend only on L or on 
other parameters that are irrelevant to our computations. Again, they may differ from line to line. 
Regarding constants K, it might appear that a constant depend on Cff 9 and C r L 9 for k ^ j. In 
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this case, we have that the constant K depends on the fc-first derivatives of and 7b and also on 
the j-first derivatives of 7 1 and 7b. We shall simplify K then into a function depending only on 
the /-first derivatives of 7 1 and 7b with / = max(fc, j), ?.e., on C^ e9 . 

From now on, we fix data ( v*,w *) and denote the associated (weak) solution w. The first step 
of the proofs is the following preliminary lemma: 

Lemma 14. Given a £ [—1,2), there exists a constant K a := K[C cvx , C^ 9 , C e u, a] s.t. vj 
satisfies: 

[ \l\ 3+a \Vw\ 2 < K a [ [ 7 «- 1 |v*| 2 +7 a - 2 | W *| 2 ] (43) 

Jb(o,l) Jb(o,l) 

Proof. Given a > —1, as 7 is a positive function, explicit computations yield that: 

A[ 7 “+ 3 ] = (a + 3)[A 7 ] 7 a+2 + (a + 3)(a + 2)|V 7 |V +1 • 

Introducing the bound (13U1) on A7, we deduce that: 

1 


|V 7 | 2 7 


2 a+1 


7 


a+2 


< 


< 


(a + 3 ) (o + 2 ) 
— A[ 7 “+ 3 ] 

^^CVX 


A [7 


q+3i 


(44) 

(45) 


Now, for any ip £ 0, L)), integrating by parts and applying Holder inequality yields: 

0< [ A[7 q+ 3 1M 2 = -2(a + 3) [ 7“ +2 V7-7>Vc 


/S(0,L) 


< (cx + 3) 


JB(0,L) 

a + 2 f 


+iF[Cr,a] /' 

Jdl 


IB(0,L) 


/ 

|^| 2 9 , 7 “ +3 da 

JdB(0,L) 


2 

+ „ 

■ [ 7 a+3 | 

a + 2 

Jb(o,l) 


2 dcr. 


ldB(0,L) 

Applying then (l44l) to bound the first term in the right-hand side of this last inequality, we obtain: 


A[ 7 a+ 3 ]M 2 


'B(O.L) 


< 


4 (a + 3 ) 
a + 2 Jb(o,l) 


[ 7 a+3 |Vy)| 2 + K[C[ e9 , a] [ I 2 dcr. (46) 

Jb(O.L) JdB(0,L) 


We multiply now (OBll with 'y a zu. This yields: 

[ 7“ +3 |VH A[ 7 “ +3 ]|zu 

Jb(o.l) ^\ a + 4 ; Jb(o.l) 


= 12 


/ B(0,L) 


w * - -div[(7 t + 7b)v*] ) 7“tn . (47) 


We compute separately the RHS and the LHS of this identity. If a < 0 the last term on the RHS 
is positive. If a £ [0, 2], we remark that vj vanishes on dB( 0, L) and apply (l46l) to w. This entails: 


LHS > 1 - 


2 max(a, 0 ) 

(« + 2 ) 


) / 7 a+3 |Vn7| 2 . 

/ Jb(o,l) 


(48) 


We note here that, since a < 2 the factor appearing in the right-hand side of this identity is 
positive. The fact that this factor changes sign when a crosses 2 makes this proof irrelevant for 
a > 2. 
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Concerning the RHS , we have first, applying (H5l) and ( 1551 ) (with <p = w and noting that 
4(a + 3)/(a + 2) < 8 ), that, for arbitrary e > 0, there holds: 


12 


/ B(0 : L) 


* c\ _ 

w 7 vj 


< 


K[C c . 


[ KlV “ 2 + |/ 7“ +3 |Vn7| 2 . 

Jb(o,l) 8 Jb(o.l) 


B(0,L) 


Similarly, by applying (l44l) . (1451) and (146)) (again with ip = vj) and (l33l) . we obtain that, for 
a £ [— 1 , 2 ) and arbitrary e > 0 , there holds: 


[ div ( 7 1 + 7 b )v* 7 a ro 

< 

a [ (7 1 + 7 b)v* • V 77 a 3 n7 

+ 

[ ( 7 1 + 7b) v * • V^7“ 

Jb(o,l) 


Jb(0,L) 


Jb(o,l) 


< 


< 


< 


K[C elu C r 2 eg ] 


K[C elu C r 2 eg ] 


K\C elh C r 2 e9 } 


I B(0,L) 


1 B(0,L) 


/ B(0,L) 


-* |2 7 a ~ 1 +e[ 7“ +1 |V 7 | 2 |w| 2 +£ f 7“ +3 |Vw| 5 

JB(0,L) Jb(Q,L) 

7“" 1 +e/ A[7“ +3 ]|ro| 2 +£ [ 7 q+3 |Vw| 2 

Jb(0,L) Jb(0,L) 

+ 7 a+3 |VH 2 - 

4 JB(0,L) 


* 12^,a—1 


This yields finally that we have, for arbitrary e > 0: 
K[C elh C r 2 eg ,C cvx \ ( f 


RHS < 


I B(0,L) 


v *| 2 7 a_1 + 


IB(0,L) 


I * |2 a—2 
\w I 7 


+ jf 7“ +3 |Vtn| 2 . (49) 

4 Jbio.l) 


/ B(0,L) 

We replace finally the right-hand side and left-hand side of ( 1571 ) with ( 155 ) 1 ) and ( 1551 ) and obtain 
the expected result by choosing £ sufficiently small depending on a. □ 

We are now in position to prove Proposition [9] and Proposition 1101 


Proof of Proposition^ Because 7 is bounded by Cq 69 on 11(0, L), it is sufficient to prove 
in the case n = 0. Under the further assumption (1341) . there holds (from the two-dimensional 
embedding H m+2 (B{ 0, L )) C C m ’ 3 ^ 4 (B( 0, L))) for arbitrary m £ N U {0}) 


|v*0r,y)| < C(|:r| 2 + |y| 2 )®||v*; U 2 (11(0, L))|| , 

\w*(x,y)\ < C'(|x | 2 + |y| 2 )5||«;*;ll 3 (ll(0, J L))||, 

for ( x,y ) S B(0,L). 

Applying then Lemma ITU to vj with a = —1/2, we obtain: 


(50) 

(51) 


f 7 S |Vtz7| 2 < K[C e i h C cvx ,C r 2 eg ] [ f 
Jb(o,l) \Jb(o, 


~s —I- W 

L) 72 72 


Here we call (l5Ul) - (l5lll and m to obtain: 

r L (rl||v*;H 2 (ll(0,T ))|| 2 , A |K; H 3 (B(0, L ))|| 2 ' 


|v *| 2 |w *| 2 


+ 


/ B(0,L) 72 


73 


5 — 


C 


ClllT 3 


(-5/2 5 
°eZZ r 


rdr 


< K[C eU \ {||v*;U 2 (H(0,L))|| 2 + |K; H 3 (H( 0 , L))|| 2 } . 


This ends the proof. 


□ 


Proof of Provosition [17/1 Because 7 is bounded by Cq £9 on 11(0,1/), it is sufficient to prove (1391) 
and (1551) in the cases n £ [0,1), n = 1 and n £ (1, 2 ). 
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Case n £ (1,2) Applying Lemma 1141 with a = n, we have, with a constant K n depending on 
C cvx ,C ell ,C r 2 e 9 and n: 


[ | 7 | 3+r, |Vw| 2 < K n 

[ 7 n_1 |v*| 2 + 7 " _2 |u;*| 2 

1 b(o,l) 

Jb(0,L) 


Here, we note that n > 1 so that 2 —n < 1. Hence, we might fixp n £ (1, oo) such that p n (2 — n) < 1. 
Introducing q n the conjugate exponent, we get: 


7 n " 2 K| 2 < 




B(0,L) \Jb(0,L) 7 (2 " )P " 

In the right-hand side of this last line, we note that standard 2D-Sobolev imbeddings yield 

||w*; L 2qn (H(0, L))|| 2 < C n \\w*; H 1 (B(0, L))|| 2 

and that, from (1311) : 

J B(0 L) 7 (2-n)p„ - K l C ell,n]j o r 2 ( 2 -„)p„ 5 #[C eH , L, »], 


as 2(2 — n)p n < 2. This ends up the proof in the last case. 


Case n = 1 Applying the embedding H 2 (B(0,L)) C C 0,3 / 4 (B(0, L)), we have: 

K(x,y)| < \w*(0)\ + C(x 2 +y 2 )i\\w*- 1 H 2 (B(0,L))\\ , V(x,y) £ B(0,L). 

The remainder of the computation follows the line of the previous case. Applying Lemma [Tdl with 
a = 1, we have, with a constant K depending on the same quantities: 


'B(O.L) 


| 7 | 4 |Vw| 2 <k 1 


,*|2 


\w 


*12 


/ B(0,L) 


< 


< 

< 


K 

I\ 

I\ 


||v*; I/ 2 (H(0, L))|| 2 + [ f H|| u ,*; jff 2 (i?( o,L))|| 2 

Jb(o,l) 7 Jb(o,l) 7 

||v*; L 2 (H(0, L ))|| 2 + £ lW h *l 0) J e ^’ + \\ w *'i H 2 (B(0, L ))|| 2 , 

IIv*; L 2 (B(0, L))|| 2 + |'ic*(0)| 2 | ln(/i)| + \\w*-H 2 (B(0,L))\\ 2 . 


Case n £ (0,1) Again, we apply LemmafHlwith a = n, and apply the embedding H 2 (B( 0, L)) C 
C 0,1 ~ n / 2 (B(O, L)). This yields with similar computations as previously: 


[ | 7 | 3+n |Vtj7| 2 < K n [ 

Jb(o,l) Je 

< K[C r 0 eg ,C ell ,n] 

< K[C r 0 e9 ,C ell ,n] 


r * |2 L ,,* |2 


B(o,L) 7 1 " 7 2 " 

*-,H 2 (B(0,L))\\ 2 rdr 


v*;H 2 (B(0,L))r 


h 1 ~ r 


"(0)| 2 r + r 3 ~ n \\w*; H 2 {B{0 1 L))|| : 


\\w*- } H 2 {B(Q,L))f 


■dr 


14 


























Case n = 0 We conclude with similar arguments, applying [Til with a = 0, and the embedding 
H 2 {B(0, L)) C C°(B(0, L)). This ends the proof. □ 


2.2 Higher order estimates: proofs of propositions [IT] and 1121 

We compute now estimates involving V fc tu with k > 1. We start with estimates away from the 
origin. 

Proof of Proposition \11\ This proposition follows from the classical regularity theory for elliptic 
problems, as developed in [10] . Indeed, fix k > 0 and 0 < e < L < oo. Applying Proposition [TOl 
with n = 2, and combining with (Kill) , we obtain for a constant K wich depends on C cvx , C e u, C^f 9 
and L that 

[Cells 2 ] 5 [ |VH7| 2 < [ |7| 5 |VH7| 2 

J B(0,L)\B(0,e) Jb(0,L) 

< K [IIv*; L 2 (B(0, T))|| 2 + \\w*; H 1 (i?(0, T))|| 2 ] . (52) 

As vj = 0 on dB( 0, L), we deduce from this inequality (applying the variant of Poincare inequality 
given in [8] Exercicse II.5.13]) that there exists a constant K := K[e, C cvx , C e u, C™ 9 , L] so that: 

\\w,H\B{Q,L) \B(0,e))[| < K [||v*; L 2 (H(0, L))|| 2 + ||«,*; ^(5(0, L))|| 2 ] . (53) 

Then, we note that the operator L := div[7 3 V] is uniformly elliptic on B(0,T) \ B(0,e) with 
ellipticity constant A £ := infa(o,L)\s(o,e) 7 3 ^ ( C e u £ 2 ) 3 depending only on C e u and e. Finally, 
applying [10] Theorem 8.8 p.183 with Theorem 8.12 p.186] if k = 0 or plD, Theorem 8.10 p.186 
with Theorem 8.13, p. 187] if k > 1 implies that there exists a constant Kff 9 depending on e, C e a 
and CIZ so that: 


|| ro ;71 fe+2 (l?(0,L)\l?(0, e ) 


< K 3 [Ik; H\B( 0, L) \ B{ 0, e))|| + ||/; H k {B{ 0, L) \ B( 0, e))||] 


< K r k e9 


l v *; H k+1 (B(0, L ))|| 2 + \\w*; H k '~\B( 0, L ))|p 


where we applied (l53l) to reach the last line. □ 

This proposition is a tool for computing the traces of derivatives of vj on dB(0,L). With this 
consideration, we obtain Proposition fl2l 

Proof of Proposition 1 121 To prepare the proof, we mention that the chain rule together with (1321) 
yields that for any integer i £ (0,..., 6} there exists a constant K[i\ depending on Cf e9 ,C e u and 
C™ 9 , for which: 

\V l ^(x,y)\<K[i\\')(x,y)f^%, V (x, y) £ B(0, L). (54) 

This inequality no longer holds when k > 6. Again, this is the reason for our lemma to hold only 
for values of the parameter k below 6. 

We prove now by induction on k for fixed m £ (0, oo) that there holds: 

” there exists a constant K^ n k depending on fc, m, C e u : C cvx , C k f 9 , with k' = max(fc + 1,2), such 
that, denoting by e*, = max(/c,3): 


'B(O.L) 


111 


3+fc+m| V fc+i ro |2 < f | 7 | 3+m |Vn 7 | 2 + \\v*-H k+1 {B(0,L))\\^ 

: L Jb(O.L) 


'B(O.L) 


\\w*-,H^(B(0,L))\\ 2 }.” (55) 
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The case k = 0 is obvious. To prove the induction argument, we fix A £ {0,..., 5} and 
introduce D k+l a (homogeneous) differential operator of order k + 1 having constant coefficients. 
Applying D k+1 to (1251) yields: 


- — div[7 3 VD fc+1 ro] = div 


R k+1 — D k+1 —( , y t + 7b) v* 


+ D k+] 


w 


(56) 


where: 


R k+1 := — (D k+1 [j 3 \7w] - 7 3 \7D k+1 w) . 


We multiply (157>1) with "f k+m+1 D k+1 oo. After integration by parts, this yields: 

— [ 1 3 + m +( k+1 '>\X7D k+1 m\ 2 = + RHS 2 , 

12 Jb(o.l) 12 


(57) 


where: 
RHS 1 = 

RHS 2 = 


f ^+ m+k d u D k+1 wD k+1 w - [ (k + m + l)V77 3+m+fc VT> fc+1 tx7£» fc+1 

JdB(0,L) JB(0,L) 

[ ( D k+1 w * - div£> fc+1 ^ + 7b)v * ^k+m+ljjk+l, 

Jb(o,l) 2 


TJ7 


VO 


f R k+1 ■ ( 1 k+m+1 \7D k+l vo + (k + m + 1)V 77 k+m D k+1 m). 
Jbio.l) 


/ B(0,L) 

+ [ R k+1 ■ u-i k+m+1 D k+1 vj 
JdB(0,L) 


We bound the first term in RHS\ by applying trace theorems and Proposition [TIJ there exists a 
constant C[D k+1 } (depending on D k+1 and also on L) for which 


\D k+1 vo\ 2 


ldB(0,L ) 


\d„D k+1 vo\ 2 ] da 


< C[D k+1 ]\\S7 k+1 vo- 0, L)/B{ 0, L/2)))|| 2 

< C[D k+1 ]\\w-H k+3 (B{ 0, L)/(B( 0, L/2)))|| 2 

< C[D k + 3 ]K r k e ^ [||v*; H k+2 (B(0, T))|| 2 + |K; i7 fc+1 (i?(0, L))|| 2 ] . 

For the second term, we apply (1321) and a Holder inequality: 


[ (k + m + l)\/^ 3+m+k \/D k+1 voD k+1 vo 

JB(0,L) 

\C e u,C r 2 9 } f ^/ 2+m+k \VD k+1 w\\D k+1 vj\ 

Jb(O.L) 

\S7D k+1 w\ 2 + K[C e u,C™ 9 } [ 7 3+m+fc| £) fc+l 

Jb(0,L) 


< K\ 


^ \ 


B(0,L) 
3+m+(fc+l) I 


B(0,L) 


To bound the last term, we note that \D k+1 vo\ < C[D k+1 ]\X/ k+1 vo\. Consequently, applying the 
induction assumption, we obtain: 

\RHSi\ < C[D k+1 }K[C elh C r 2 e9 ,l CT,^+i]{ / |7| 3+m |Vtn| 2 + ||v*; H k+2 (B(0 7 L))\\ 2 

k Jb(o,l) 


+ \\w*;H ek + 1 (B(0,L))\\ 2 }+ ^ f 7 3+ra+ ( fe+1 )| V £) fc+1 ro| 2 . 

’ 2 J B (o,l) 
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As for RHS‘ 2 , we distinguish between two cases. If k + to > 1 we have that 2(m + k + 1) > 
3 + m + k so that a standard Holder inequality yields, with the induction assumption: 


[ ( D k+1 w * - idivD fc+1 ( 7t + 7 {,)v*)7 fc+m+1 D fc+1 n7 

JB(0,L) 2 


< C[D k+1 } 


||w*; H k+1 (B(0 1 A)) || 2 + CIS ll v *; H k+2 (B(0, L))|| 


< C{D k ^]K{K s ™°,C r k e ° 2 ] 


[ 7 3 +m|y a7 |2 _|_ || v *; H k+2 (B(0, L))\\ 

Jb(o,l) 


B(0,L) 
2 


^3+m+fc|^7A;+l„|2 


W\ 


||w*; H ek (B(0 1 T))|| 2 


Whereas, if 0 < [k + m ) < 1 so that k = 0 in particular, we introduce arbitrarily the missing 
powers of 7 and recall (l32|) - (l33)) : 


[ ( D k+1 w * - idivD fc+1 ( 7t + 7{,)v*) 7 fe+m+1 D fe+1 n7 

Jfl(o.i) 2 

|Vu>*| 


C 


1 (fc + m) 

IB(0,L) \ 'y 2 2 


7 2 2 i=o 


( fc + m ) I 3 JL 1 1 

, 2 +2 D + W 


where 1 — (fc + to) < 1 so that: 


's(o,i) 7 


< ^[C' en ]||u>*; J H 3 (B(0,L))|| 2 , 


and 


/ 

. L., Edivvi 

Jb(o,l) 

7 2 2 j =0 


< K[C e i h C r 0 e9 }\\v*-H 2 (B(0,L))\\ 


Hence a Holder inequality yields that: 

[ ( D k+1 w * - idiv£> fc+1 ( 7t + 7b )v*) 7 k+m+1 D k+1 w 

Jb(o,l) 2 

7 3+m |Vn7| 2 + ||v*; H k+2 (B(0, L))\\ 2 + ||to*; H ek (B{ 0, L))|| s 


B(0,L) 

<K[K^,C r 2 ea ,C ell } 


/S(0,L) 


To treat the last terms, we expand the differential operator R k+1 and apply (l54l) yielding that 
there exists also a constant C[D k+1 ] for which: 


|i? fc+1 | < C[D k+1 ]Y^ |V i+1 -ccr||V fc+1_Z 7 3 | < C[D k+1 ]^TK[k+ 1 - Z]| 7 | s 


— |V‘ + i T77|. 


1 -0 


1 -0 


Consequently, we bound similarly as above by applying trace theorems and Proposition IIII 


f R k+1 ■ v^ k+m+l D k+1 vo 
JdB(0,L) 


< C{D k+1 }K[C r k e ^C elu C r 2 e9 ]\\w- H k+1 (dB( 0, L))|| 2 


< C[D k + 1 ]K[C r k '* v C M , C r 2 eg \ ||t*7; H k + 2 (B( 0, L) \ B{ 0, L/2 ))|| 2 

< ClD^KiC^Ceu,^ 9 ] [IIv*; H k+1 (B(0 1 A))|| 2 + ||w*; H ek (B(0 1 A))|| 2 ] . 
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For the other term, we recall (1521) and apply the induction assumption. This yields: 


[ R k+1 ■ ( 1 k+m+1 \7D k+1 w + (k + m + l)\/^ k+m D k+1 w) 

J B(Q,L) 


< K[C r k l\,C elh U 2 -’\2^ I 7 

1=0 Jb(o,l) 


■|v i+1 n7| ( 


3+m+(fc + l) , . . 3+m+fe . i i A 

7 2 \X7D k+1 m\ + 7 2 |V fc+1 iu|J 


< 




'B(O.L) 


7 


3+m+(fc+l )\\J 


+K[C% X ,C M ,C?°]K™° 

+\\w*;H^(B(0,L))\\ 2 


' B(0,L) 


| 7 | 3 + m | Vttj| 2 + ||v*; H k+1 (B(0, T))|| 5 


We obtain finally that there exists a constant K^ n £ +1 depending on Kj ..^, K s ™'£, K^ 9 and 
C r k e + % Cr,C ell for which: 


RHS i 


+ RHS 2 


12 

< C[D k+1 ]KZ n k g +1 


' B(0,L) 


| 7 | 3+m |V W | 2 + ||v*; H k+2 (B(0, T))|| 2 + |K; H*™ (B(0, T))|| 2 


+ l. 


B(O.L) 


7 


3-|-m-f-(/c+l) 


|V£> fc+1 T 


We introduce this inequality into (f57l) and obtain that: 

[ |^|3+m+(fc+l) |y £) k+1 w\ 2 

JB(O.L) 

| 7 | 3+m |Vw| 2 + || v* ; H k+2 (B(0, L))\\ 2 + \\w*; H ek+1 (5(0, T))|| 2 


'B(O.L) 

< C[D k+1 ]K™° + 1 


lB(0,L) 


Given the dependencies of and . the operator D k+1 being arbitrary, we obtain 

the expected inequality for the rank k + 1. This ends the proof. 

□ 

2.3 The case of spheres. 

We end this section by computing sharp asymptotic expansions of 

[ 7 3 |V W | 2 

Jb(o,l) 

when 7 represents the distance function between one sphere of radius S and another one of radius 
R. Until the end of this section, we assume that, close to the origin, we have: 


7 y) = S - ^/S 2 - (x 2 + y 2 ), j b (x , y) = -R + \JR 2 - ( x 2 + y 2 ). 
Hence, 7 depends on r = \Jx 2 + y 2 only and we have the Taylor expansion: 


2 


=-K O(r^), ; ; 

2s r 2 ’ l{r) = h + i^- + i^3 +°( r6 ), Close to r = 0, (58) 

7 M = -Ijj+oe 4 ), 2Ri 8B » 


18 



















where Ri and R 3 satisfy ( 1 ^ 1 ) . 

Then, we split f = fo + fi + Ir, where: 

fo = w*{ 0 ), f 1 = 1 -(^-^Jv*( 0 )-e r +xd x w*( 0 )+yd v w*( 0 ), f R = f - (/ 0 + / 1 ). 

Due to the linearity of the divergence problem (I26p - (I27[) . the solution w admits a corresponding 
decomposition: w = wq + wi + wr with obvious notations. In this section, we compute separately 
the asymptotic expansions of the quantities 

[ 7 3 |Vtu 0 | 2 , f 7 3 |V^i| 2 , [ 7 3 |VtuR| 2 . 

Jb(o,l) Jb(o,l) Jb(o,l) 

First, as f R corresponds to / from which the first orders in the Taylor expansion around 0 are 
subtracted, we may introduce £ C°°(B(0, L); R x R 2 ) s.t.: 

fR = w* R ~ ^div(( 7t + 7 b )v * R ), 


We have then: 

«4(O) = 0, V*, y ti$(0) = 0, vjj(0) = 0. 

Hence, Proposition [9] entails that, for h £ (0,1] : 

/ 7 3 | Vn7f {| 2 = 0(1) {\\u*;H 2 (B(0,L))\\ 2 + \\w*-H 3 (B(0, L))\\ 2 } . 

J B(0,L) 

Here and in what follows, we denote by 0(1) a quantity which depends on (R , S, L ) and h and is 
bounded by a constant depending only on R , S , L whatever the value of h £ (0,1]. 

It remains to treat the two cases of vjq and w\. We remark at once that, for wq, the associ¬ 
ated source term fo is constant, while, for zuj , the associated source term reads fi = fi(r,6) = 
f c r cos (9) + f s r sin( 0 ) where: 

fc = \ (ij “ ■ Gl + f s = \ v *(°) ■ e2 + dyW*{0). 


We start with wq: 

Proposition 15. Under the assumption that 7 is radial and satisfies 

R\ 3Rf 


/B(0,L) 


7 3 |V G 7 0 r =72 tt| / 0 r 


Rl 


IM»| 


, there holds: 

+ 0 ( 1 ). 


(59) 


Proof. By linearity, we only treat the case fo = 1. Under our symmetry assumptions, the unique 
solution wo to (I26f) - (127f) with / = 1 is certainly a radial function and explicit computations yield 
that: 

6 s 6r 

^)ds, drwo(r) =-^y, Vr £ (0,L). 

Consequently, we have: 

9 f L r 3 dr 

T 3 Vu7o 2 = 72n / 

,L) Jo 7 3 U) 

At this point, we note that 

r 2 r 4 

70) =h+ 2R 1 + 8 R 3 + rem ^ 
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with |rem(r)| < Cr e for all r £ (0 ,L). Consequently, introducing r o small enough so that we 
might expand 7 in power series, we have: 


l 


L r 3 


7 ( r ); 


i-dr = 


7 (r ) 3 
r ° r 3 dr 


dr + 0(1) 


r r 0 


r' dr 


/0 (^ + ^r) 3 8 R 3 J° ( h + ^ 4 


0 ( 1 ) 


1 f Vh s 3 ds 


3 f Vh s 7 ds 


(i + ifr ) 3 8 S U° (1 + *) 


0 ( 1 ), 


2R X ■ 


where explicit computations yield: 

rJk s 3 d s 


(i + ^) 3 


= R{ + hO{ 1) 


■Jk s 7 ds 


(i + ^r) 4 


= 8 i?f | ln(/i)| + 0 ( 1 ). 


Finally, we obtain: 


[ 7 3 |Vn7 0 | 2 = 72it 

J b(o,l) 


R\ 3 R\ 
~h ~ “R 3 " 1 


ln(/i)| 


+ 0 ( 1 ). 


□ 


We end this section with the case of an /i-like source term: 

Proposition 16. Assume that 7 is radial and satisfies (1581) . Given a source term of the form 
f{r,0) = f c rcos(9) + / S rsin(0) on B(0,L), there holds: 

f 7 3 |Vnf = (|/ c | 2 + |/ s | 2 ) ln(/i)| + 0 ( 1 ), (60) 

Proof. Up to shift 6 with a phase and call the linearity of our problem, we prove the above result 
in the case f c = 1 and f s = 0. Then, the proof is divided into 3 steps. 


Step 1: reduction to an ode. Under our symmetry assumptions, the unique solution to (l26l) - 
(l27l) reads w(r, 6) = qn{r) cos {6) where qh is the unique solution to the ode: 

^ d r [ 7 {r) 3 rd r q h ] - 7 - —12r, r £ (0, L) (61) 

q h (r) = 0, r £ {0, L} . (62) 

The boundary condition in r = L is the translation of w{r, 6) = 0 on r = L while the boundary 
condition in r = 0 is derived from the condition that: 


7 r 



\d r qh(r)\ 2 


\Qh(r) I 2 


rdr 


' B( 0 ,L) 


7 3 |Vn7| 2 da;dj/ < 00 . 


We note that (ill) rewrites: 




















Step 2: Construction of an approximate solution. We introduce then q the unique solution 
to an auxiliary problem. The construction and properties of this auxiliary function are stated in 
the following lemma whose proof is postponed to the appendix: 


Lemma 17. Given R\ > 0, there exists a unique q £ C([0,oo)) D C'°°((0,oo)) solution to 



g(0) = 0 lim q(s) = 0 . 

s—yoo 


Furthermore we have the asymptotic description: 


48 Rl 144 Rl rem(s) 

q(s) = + rem(s) d s q(s) = — 


5s 3 

where |rem(s)| < K[Ri]/s a for s > 1. 

With that auxiliary function at-hand, we set: 

L 


5s 4 


Qh O') = tt? ( Xr) “ TX7T9 


, Vr€(0,L), uj = qh{r) cos(9) on 13(0, L). 


hVWhJ tfhVWh; 

Introducing y(r) = (h + r 2 /(2??i)), we obtain by substitution that qn is a solution to 

^ d r [7 {r) 3 rd r q h ] - ^ = -12 r - jjjtQ x( r ) , rg(0, L) 

q h (r)=0, r £ {0, L} . 


(63) 


where 

X{r) = ^ d r [y(r) 3 r 2 ] -7 3 (r). 

Consequently, lo is an ^-solution to: 

-^ div [7 3 Vro] = (r + g X( r )^ cos ( e ) > on L )> 

w = 0 on r = L . 


We note that, for n > 3, there holds: 

L := 


f 7 " Vot 2 da;d y = 7r 

f L l( r ) n 

10 - ( M2 . \Qh{r)\ 2 ' 
\drqh{r)\ + 2 

rdr 

1B(0,L) 

Jo 

r 



Replacing q\ with its value and noticing that ^ = 0(1) because of (l63l) . this entails: 


L 2 m 


In = 


h 3 

+0 


\ru + ^y 

\d r q(r/Vh )\ 2 \q(r/y/h )\ 2 

rdr 

Jo V / 

h r 2 



1 + 


2R 1 


\d r q{r/Vh )\ 2 \q{r/Vh )\ 2 

j I O 


rdr 


We bound the first integral on the right-hand side by changing variable r —> \fhs. This yields: 
| d r q{r/Vh )\ 2 \q(r/Vh )\ 2 


rL ' r 2 V 


h + 


2RiJ 


rdr 

"L/y/h 


= h n 


1 + 


21?! 


\9 s q(s)\ 


2 , l^)! 2 


sds. 
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Given the asymptotic expansion of q this entails: 


h + 


2 i?i 


\d r q{r / y/h)\ 2 \q{r/y/h){* 


(■ r/Vh ) 2 


h 3 R ? 


rdr = < 40 

h 3 Q( 1 ) 


1 (48 2 |ln(/i)|+ 0(1)) if n = 3 
if n > 4. 


In particular we obtain finally that: 


288R 3 i 


[ 7 n |Vro | 2 = < 5 

JjB(0,L) I 0(1) 


| ln(/i)| + 0(1) for n = 3 

for arbitrary n > 4. 


(64) 


Step 3: Computing the distance between vj and vo We introduce the difference D[w\ := 
vj — w. It is the solution to: 


- Idiv^VDM] = -4div [+ - 7 3 )Vro] - Ii + J , (A 
D[oj\ = 0 

Multypling ([65jl with D[oo\/') 1 we obtain, after integration by parts: 


x(r)cos(0), on 6(0,6) (65) 
on 9-6(0, 6 ) ( 66 ) 


1 

12 


/ 7 2 |VO[tu]| 2 -i [ V| 7 | 2 V|OH| 

Jb(o,l) 4 Jb(o,l) 


>B(0,L) 

On the left-hand side of this identity, we have: 


= RHS. 


(67) 


1 

12 


/ 7 2| V OH| 2 -i/ V|7| 2 V|OH| 

Jb(o,l) 4 Jb(o,l) 


B(0,L) 


1 

12 

1 

> — 

“ 12 


1 


/ 7 2 |vohi 2 + -/ [A| 7 | 2 ]|oni 2 

Ib(o,l) 4 Jb(o,l) 


'B(O.L) 


7 2 |VO[td7]| 2 


/S(0,L) 


7 | 0 [ n 7]| 2 


as A| 7| 2 = 2[yAy + |Vy| 2 ] > 2cy on 6(0, 6 ). On the right-hand side, we obtain the bound above: 


\RHS\ < 4 + 


12 126 2 /ii 


where, after integration by parts: 


h = 


< 


f div [(y 3 — y 3 )Vti 7 ] 
Jb(0,L) 


D[z 


lB(0,L) ~ ~ 7 

( |73 ~ 73| 1V7||V^||6M| + 


/ B(0,L) 7 


r | 7 3 -7 
B(0,L) 7 


- 1 V tzj 11 V Z7 [txj] | 


C 

< — 


(7 3 -7 3 ) 2 /l V 7 | 2 , 1 


e 7 b ( 0 ,l) 7 2 V l7l 3 


+ A) l v ^| 2 + £ / (y 2 |V 6 >[w ]| 2 + 7 | 6 [tu]| 2 ). 

7 / Jb(o,l) 


for arbitrary £ > 0. At this point, we remark that, as 7 and 7 share the same (non-vanishing) 
Taylor expansion up to the second order, we have, with a constant K depending only on R, S : 


I 7 3 — y 3 | < A'y 4 . 


22 
























Recalling ( 1521 ) to bound V7 and ( 1551 ) we obtain: 

f (7 --77 nvzj! iw. 


Jb(o,l) 7 2 V 111 3 7 

and finally, for arbitrary positive e, we have: 


4-^)1 Vw\' Z <Kl 7 4 |Vw| 2 = 0(1), 


'B(0,L) 


r 0 ( 1 ) 

h < —j 1 + £ 


B(0,L) 


(7 2 |VD[ro]| 2 +7|D[n7]| 2 ). 


Similarly, we bound 4 : 

h = 


X{r) cos(6>)^^i 
B(0,L) 7 


C f L 

5 jJ. 


-9,. [7(r) 3 r 2 ] - 7 3 (r) 


rdr 


+ £ 


7WI 2 , 


7 J B(a,L) 

where there is a constant AT depending only on the characteristics of 7 for which 


Hence, we have: 


-d r [7(r) 3 r 2 ] - 7 3 (r) 


-d r [7(r) 3 r 2 ] ~7 3 (r) 


< ky 


rdr 

— = 0 ( 1 ) 

/yO 


and finally, for arbitrary s > 0, there holds: 


2 , 0(1) 


h<£ l\ D [w}\ 

JB(0,L) e 

Introducing the above computations of the left-hand side and right-hand side into (1671) . noticing 
that 


1 


12 L 2 hi Wh 
thanks to (l63l) . we obtain, choosing e sufficiently small: 


= 0 ( 1 ) 


[ 7 2 |V£M| 2 +/ l\D[w }\ 2 = 0(1). 

J B(0,L) Jb(O.L) 


B(0,L) 


( 68 ) 


Step 4 : Conclusion. We are now in position to compute the asymptotics of 


4 := / 7 3 |Vw| 2 . 

Jb(o,l) 


Indeed, there holds: 


4 = 


f 7 3 |Vw| 2 + f [ 7 3 -T 3 ]|Vtz7| 2 
Jb(o,l) Jb(o,l) 

+ [ 7 3 |VL>[ro]| 2 + 2 f 7 3 VD[ro] : Vro 

Jb(0,L) Jb(0,L) 

= [ 7 3 |Vw| 2 + Ji+ J2 + J3. 

JB(0,L) 

As previously, we bound I7 3 — 7 3 | < ATy 4 so that recalling (l6dl) with n = 4 we have: 

|*4| < K f 7 4 |Vd7| 2 = 0(l). 
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Then, (1551) implies: 


\M <K f 7 2 |VD[n7]| 2 = 0(1). 

Jb(o,l) 

We conclude with similar arguments and applying Proposition [TO] in case n = 1 (here vj is a 
solution to flUMZZD with a special v* and w* =0): 


\M = 


< 


I y 3 V D[vj\ : V vo 
Jb(0,L) 

[ 1 2 \VD[w}\ 2 + [ 7 4 |Vw| 2 = 0(1). 

Jb(o,l) Jbio.l) 


Hence, there holds: 


h = 


[ 7 3 |Vd7| 2 + 0(1) = ln(/i)| + 0(1), 

J B(0,L) u 


because of (1M1) . This ends the proof of Proposition [TCI 


□ 


3 Proof of Theorem |4] 

In this section, we complete the proof of Theorem [4] We first introduce some notations: 

• \l is a truncation function which vanishes outside B( 0, L). Namely, we set 

XL(x,y) = x > V (x, y) £ R 2 , 

with x G C(?°(R; [0,1]) satisfying 

X (t) = 1, Vt € [-1/2,1/2], X (t) = 0 , Vt G R \ [-1,1]. 

• we recall, that given £ > 0, we denote: 

Ge = {(x,y,z) G M 3 s.t . \(x,y)\ < £ and z G (yb^’, y), h + 7 t (a:, y))} , 

• Q* is a smooth subdomain of \ Q L / 4 which does not depend on h and satisfies: 

fP C T whatever h £ (0,1], dfl \ 8 Gl/i C dfl*. 

We note that such an f2* exists as, thanks to assumptions (A2) and (A3), there exists 8 > 0 
depending on 6 and C e ii such that, whatever the value of h G (0,1] there holds 

{(*, y,z) s.t. dist((a;, y, z ), dSl) <S}\ Q L/4 C T. 


3.1 Construction of asymptotic approximation 

Let fix one boundary condition u* £ C°°{dGl) satisfying (l2lT) . At this step, we explain how to 
construct a velocity-field v[u*] which shall approximate the solution to the Stokes problem with 
boundary conditions u* £ C°°{dGi) in the regime h « 1 . 

Let split u* into its tangential and normal parts (according to the tangent space in the origin): 

u* = u* + u* ± e z £ C°°{dn). 
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We introduce qi n the unique solution to 

-pj div (7 3 V<7in) = u* ± - ^div( 7 t + %)u* - ^(h - 27 6 )divu* , on 5(0, 5), ^ 

qin = 0 , on 95(0, L). 

We recall that the data 7 t, 7 f, and u* are smooth so that we have existence of a unique qi n solution 
to this equation satisfying: 

qin £ C°°(5(0, 5)). (70) 

Second, we construct an auxiliary pair velocity-field/pressure in the aperture domain. We denote: 

Pin{x 7 y) = XL(x,y)q in {x,y ), V (x,y, z) £ Gl . (71) 


and 


Vin,//(x,y,z) 


Vin,±(x,y,z) 


-(z- {h + ^t)){z ~lb)^x,yPin + 


h + 'yt- z 


XL Uy 


^div x>y 


rh+j t 

/ (s-(h + -y t )){s-jb)^x, y Pinds 

J Z 


(72) 

(73) 


+ 




Xl-u-h 


ds, 


for ( x,y,z ) £ Gl- We note that this velocity-field vanishes outside Gl■ We keep notations to 
denote its trivial extension to T in what follows. 

Combining (1701) with (R71 - (1771) - (1731) . we obtain that: 

(VR1) \in £ C'°°(J r ) and has support in Ql, 

(VR2) div x in = 0 on J 7 , 

(VR3) x in = 0 on dB , 

(VR4) xi n = u* on dGl n dQ L / 2 

All properties satisfied by x. m are obvious but : 

Vin,x{x,y, 7 b{x,y)) = u* ± {x,y ), V (x,y) £ 5(0,5/2). 


This property is a consequence of qi n solution to (1M1) . We warn the reader that, to check this 
implication, one cannot commute the divergence and integral operators prior to taking the trace 
onz = 7 b(a :,y). Indeed, the differential operator div X)y does not include only derivatives parallel 
to the boundary z = 7 b(x,?/) so that one has to expand the differential operator div X) y inside the 
integral to perform this computation. 

The internal velocity-field Vj„ does not match the boundary condition u* outside Gl /2 (unless 
u* vanishes). As we expect no singular behavior of the exact solution to the Stokes problem far 
from the singularity, we extend the approximation in a simple way outside Gl/ 2 - Precisely, we 
define 

x* ext = u* - Vj n |an on dGl. 

As x in £ C°°(F) and u* is smooth, we have that x* xt £ C°°{dGi). Moreover, the properties of 
the boundary values of Xi n that we mentioned above ensure that x* xt vanishes on dGl D OGl/ 2 - 
Consequently, we might extend x* xt by 0 to dGl*, defining in this way a smooth function. Then, 
we denote ( x ext ,q ex t ) the unique solution to: 


ext ^ Qext 

= 0 , 

on Gl* , 

(74) 

V • V ex t 

= 0 , 

on Gl* , 

(75) 

V ext 

* 

V ext 5 

on dGl* . 

(76) 
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This solution indeed exists as: 




• ndcx 



• nder — 



• ndu 



nder = 0. 


As dfl* and v* xt are smooth we also have that v ext is smooth on and vanishes on <9f2* n9[J r \n*]. 
Hence, the trivial extension of v ex t to T (that we still denote v ex t for simplicity) satisfies: 

(VR5) v ext is continuous piecewise smooth and has support in Q* 

(VR6) div v ext = 0 on T 


(VR7) v ext = 0 on dB, 

(VR8) v ex i = u* - v in on dfl 
Our asymptotic approximation v[u*] reads then 


v[u*] = v in +v ext . 


(77) 


Even though the vector-field v ext is important in order to obtain a velocity-field v[u*] that 
matches the boundary conditions u* on the whole dCl, this external velocity-field does not contain 
any information on the way the velocity-field v[u*] diverges when h is small. Indeed, we have the 
following lemma: 


Lemma 18. 

such that: 


There exists a constant K depending on C™ 9 ,C e u,dfl and 5 (introduced in (A3)) 


|Vv R 


< K 


IT 


|u*; iJ 2 (i?(0, L))|| + ||u*; iJ 2 (<9fi)|| 


i 2 


Proof. Indeed, as solution to the Stokes problem on H*, the exterior velocity-field v ext satisfies 


IT 


|Vv ea , t | 2 


= / |Vv ext | 2 

Jn* 

< C*^w* ext \(dVL*)\\ 2 . 


where the constant C* depends only on Q* which depends itself on the geometry of dQ away 
from the singularity and S. Moreover, recalling that v* xt vanishes in Q L / 2 , we might construct 
a function f £ C 00 (H*) such that ( = 1 on Supp(v* xt ) and ( = 0 on cAT \ dfl. We have then, 
introducing T[u*] £ H 1 (fl) a lifting of u* on the whole Q: 

||v^ t ;JTi ( an*)|| < ||Cu*;fri(afi*)|| + ||Cv in ;^(an*)|| 

< ClCr[u*]; J H 1 (H*)|| + ||Cv m ;i7 1 (H*)|| 

< C f *(||u*;^(SH)|| + ||v i „;iJ 1 (fI*)||) . 


As Q* remains away from the singularity, we might bound from below 7 by a constant 6 depending 
only on C e u on Ql D O* and apply the explicit formulas for Vi n to obtain the following bounds: 


||v m ; H\n*) || < C*(l + 1 ) (||7&; C\B{ 0, L))|| + || 7t ; C 2 {B{ 0, L ))||) x ... 

... x (\\q in ;H 3 (B(0,L)\B(0,L/2))\\ + \\u*;H 2 (B(0,L))\\). 
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We now wish to apply Proposition [TT] to qi n in order to yield a constant that bounds its U 3 -norm 
by something which does not depend on h. However, one term in the right-hand side of (1691) 
depends on h. To avoid this difficulty we expand qi n = q^ + hq' 2 J with pressure field components 
defined as respective solutions to: 


^-div(7 3 Vg£ } ) = 01 + 7&divu / * ) - idiv( 7t + 7 b )u* , on B{ 0, L ), 

Qin = 0, on <95(0, L). 


(78) 


and 


— J_div( 7 3 V^ } ) = -idiv U/ *, on 5(0, L ), 


(2) 

Qin 


= 0 , 


on dB( 0, L). 


(79) 


In particular, we remark that q^ is a solution to with a source term / given by (1281) 

associated with 

w ± = u ± + 7 &divu / * v* = u* , 


and that 


is a solution to ESl) - (l271) with a source term / given by C51) associated with 


w]_ = — -divu 7 * v* = 0. 

Hence, we apply Proposition [TT] to both q^J and q^ and obtain: 

\\q%; H 3 (B(0 7 L) \ 5(0, L/2))|| + h\\q%>; H 3 (B( 0, L) \ B( 0, L/2))|| < Ju*; H 2 (B(0, L ))||. 


This entails finally: 


||v m ; H\Q *)|| < C* (l + 1) (||7fc; C 2 (5(0, T))|| + || 7t ; C 2 (5(0, i))||) x ... 

...xK[ eg ||u*;7T 2 (5(0,5))||. (80) 

□ 


3.2 Main steps of the proof of Theorem [4] 

We fix now u* £ C°°(dfl) satisfying (I2T1) as in the assumptions of Theorem [4] We split this 
boundary condition in 

u * = U L + U *R 

with: 

U L = u*±(0)e z + u* (0) + (xd x u*j_(0) + ydyU*j_{0))e z , 

u* R = u* - u* s . 

Note that, straightforward computations entail: 



• ndcr = 0, 


u*er(ffl), Vie (as, 5}. 


So, we might define the solutions of Stokes system Jl]) ([4|(-(-([6])) as given by Theorem [3] for 
boundary data u*,u* s ,u]j. We denote these solutions (u ,p), (u as ,p as ), (u r,Pr) respectively. We 
also introduce the asymptotic approximations: 


v as = v[u* s ] v H = v[u;y. 
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Let denote v = v as . Then, as u* s ^ is constant, we have that 

(h - 27 b)divu* Sj// = 0 

so that the pressure associated with v os satisfies m-m and, consequently, v as satisfies m-m 
in Sl/ 2 ■ To complete the proof, it remains to compute 

IIu v; V\\ = ||u aa + - v QS ; V\\ < ||v as - u as ; K|| + ||u fl ; V\\. 

Next subsection is devoted to the computation of || u^; V|| and the following one to ||v as — u as ; V\\. 
In particular, the proof of Theorem [I] is completed by applying Proposition [TD] and Proposition 
I2T1 in the general case and by applying Proposition [191 Proposition [22] and Proposition [23] in the 
radial case. 


3.3 Asymptotics of ur 

We start with the remainder term u#. Keeping notation wr for v[ujj], we prove: 

Proposition 19. There exists a constant K depending on <912, <5, C^ eg , C cvx , C e u and L such 
that, if h G (0,1], there holds: 


it 


|Vufl| 2 < K \\u*;H 3 (B(0,L))\\ 2 + ||u*; ffi (cK2)|| 5 


We recall that thanks to the variational characterization Proposition [G] u/v realizes the mini¬ 
mum of V norms among divergence-free velocity-fields w satisfying the same boundary conditions 
as u r (i.e. w G K[uJj] with our notations). By construction, we have wr G K[u|j] so that the 
above proposition is a consequence of the lemma: 


Lemma 20. There exists a constant K depending on on dfl, 5, C ^ eg , C cvx , C e u and L for which, 
if h£ (0,1] : 


J |Vv fl | 2 < K [||u*; H 3 (B(0, L ))|| 2 + ||u*; ff*(0n)|| s 


(81) 


Proof. We drop index R in v for the whole proof. We recall that, by construction, we have 
v = Wi n + vext where v,„ and w ext are computed via d-d and d-d respectively. We first 
apply Lemma [lH] showing that we only have to focus on the following contribution of w ln : 


/ |Vv in | 2 = / |Vv m | 2 

>T JQl 


Explicit computations show that: 


|^ x,yVin, II 


— Cl ''f ' x,yPin\ T" 7 \^x,yPin I + 

i=t,b ' 

\dz^in,//\ < C L (^\\ 7 X} yP in \ + ^-\u* R// \j , 


| Vx,y7i 


+ l u k//l + \^x,y n *Rj, 


— Cl x,yVin,// | ? 

and: 

I Va:,yUm,_L | < C L \ ((I '^ x,y7i + I' ,yH\j 2 ) I '^ x,yPin \ + l 2 1' Vx,y7i \ \^l, y Pin |) 

V i=t,b 

+ ( |V 2 , y7t | + |V 2 >y7 | + |V x ,j/7t| + \Vx,yl\ + l Va! -»Ttll v »,«7l + |Vx,y7l 2 + 7 . ]U * R // 1 


+ (|Vx,„7t| + l V 7l + 7) \Vx,yU* Rj/ I + j\V 2 XyV u* R // 1 + T\V 3 Xty p in \ y , 


28 








with Cl a constant depending on L. Introducing that, for all (x, y) £ 5(0, L ), there holds 

I'V*,„76| + + |Vx lW 7l < C M }^ , |V^ y7i | + |V^ 7 | + |V^ 7 | < K[C r 2 ea ], 

we obtain that, on Ql, there holds: 

|Vv m | < K (j\W x ,yPi„\ +'y 2 \'Vl, y Pin\ +l 3 \Vl,yPin\ + ^1 + ^ |u^ f// | + |'H I + 7 |' Vl t yU* R> , 
Integrating this inequality entails that: 


>Gl 


\X7v m \ 2 dxdydz < K[C ell ,C r 2 e9 ,L}\\u* R// ;H 2 (B(0 1 L))\\^ 


+ K[C ell ,C r 2 e9 ,L\j^^ ^7 3 |V* >y p in | 2 +l 5 Wl y Pin\ 2 + 7 7 |V^p iri | 2 + ^^j dxdy. (82) 

In this last identity, we note that, for k £ {1, 2, 3} : 

k -1 

\^t, v Pin{x,y)\ < \^ k Xty qin(x,y)\+C L '^2l Bi0:L) \ B{0tL/2 ){x,y)\W J xy q in {x,y)\, V{x,y) £5(0,5). 

j=o 

So, we split again qi n = q^ + hq^ with q^ and q^ defined respectively as the solutions to 

~1 ~ 2 dWx ’y(' y3 'Vx,y ( lin) = ( U R ,-L + 76 div Xiy u^ )/r/ ) - ^div Xjy ( 7t + 76 ) u* R // , on 5(0, L ), 

on dB{ 0, L). 

(83) 


and 


Qin = 0 . 


1 ’• ' t (2)n _ 1 


(84) 


-^ div ^(^ 3V *.^) = -^ div a:, 3 / U fl,// . on 5(0, L), 
qin = 0, on 95(0, L). 

So we have that q^ is a solution to (l26l) - (l27l) with a source term / given by (l28l) associated with 


W ± = -~ A ™ X ,y\l* Rj/ 


v* =0. 


where we note that ||div Xiy uJj ^; L°°(5(0, i))[| < C||uJj; 5 3 (5(0, L))[|. Applying propositions HT1 
[Hand[TO] to u- 2 \ we obtain that, 


[ {l^x^qfnl 2 + 7 5 |V 2 , !/ 9in ) | 2 + 7 7 |V Xi!/ ^ ) | 2 ) dxdy 

J B(Q,L) v ' 


+ ||^; 5 3 (5(0, L) \ 5(0, L/2))[| < 


A1|u^;5 3 (5(0,L))||^ 


with K depending on C ^ 9 , C c „ x , and 5. 

As for q\n > we remark that it is a solution to m-m with a source term / given by 
associated with 

w ± = u *r,x + 7&div x> 2 / Ufl i// v* = n* Rj/ . 

In this case, we have that ud[_(0) = 0 and also that ViAj_(0), v*(0) vanish. Consequently, we are 
in position to apply propositions ITT! [TO1 and [51 to q^, yielding that: 


/ (7 3 |V,. y <,«| ! +7 s |V^ 5 «| 2 + 7’|V^ 9 0»| 2 ) 

+ \\q£>; 5 3 (5(0, L) \ 5(0, L/2))|| < tf||u£; 5 3 (5(0, L))|| 2 , 
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with K depending again on Cff 9 , C cvx , C e u and L. Combining the computations for and , 
and arguing that 

IIujj; H 3 (B(0, Z/))|| < C\\u*; H 3 (B(0, L))|| 


we get hnally: 



with K depending on C ^ eg , C cvx , C e a and L. 

For the last term on the right-hand side of (1571) . we add that u,* (0) = 0, implying: 


I u *R,//( x ,y)\ < (M + M)||u* ;j ff 3 (i?(0,L))|| , on B(0,L). 

Consequently, going to polar coordinates yields (recall 7 satisfies (13T1) whatever the value of h £ 

( 0 , 1 ]): 


R,//\ K 27r||u*; U 3 (.B(0,1/))|| 5 


IB(0,L) 7 


Cell 


r Z Cell 


This ends the proof. 


□ 


3.4 Asymptotics of v as — u as 

We proceed with a first bound on the singular term in the general case (z.e., without structure 
assumption on 7 ). We prove: 

Proposition 21. If h £ (0,1], there exists a constant K depending on 90, 5 , C^ eg , C cvxj C e u and 
L such that there holds: 


J |Vu QS - Vv as | 2 < A'[K(0)| 2 |ln(h)| + ||u*; tt 3 (H(0,L))|| 2 ]. 

Proof. We recall that v as = Vi n + v ext . It is sufficient to compute a constant K (independent of 
h) such that, if h £ (0,1], for all w £ Vo, we have: 


(A w in - v Pin) ■ w 


<T 


< K 


|Vw|^ 


it 


(with pi n given by (1711) 1. Indeed, we have then: 


|(Av qs - S7pim w)| < 


J (Avj n - \7p in ) • w + J ' 


Vv„,( : Vw 


Applying Lemma [18] together with the remark that the mapping u* 1 —> u* s is linear continuous 
H 3 (B( 0, L)) '? H 2 (B( 0, L)) (~l we bound the other term by: 


Vv„ f : Vw 


IT 


<^||W ;j ff 3 (I?(0,L))||||w;V||. 


with K depending on 90, C^ 9 , C e u and S. 

We emphasize that v,; n and pt n have support in Ql so that: 


I[w] := / (Av m - Vp in ) • w = / (Av in - Vp in ) 
JT JQl 


By construction, we have that: 


^zz^inj/ — V x,yPin , 9 z Pin — 0 
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Consequently, there holds: 


^ Pin — ^x,yVin,// “ 1 * ^^in.-L^z * 

For any w € Vo we can then bound 7[w] by integrating by parts: 


|/[w]| = 


/ ^x.yVin,// :V x , y W + / Vu in ,J. • Vffll 
IQl JQl 


< 


' Gl 


\^x : yVin,// j F|VUi ni J_| 


|Vw;L 2 (J-)||. 


The remainder of the proof follows the line of the proof of Lemma EDI First, we bound V ; , : , y v in;// 1 
and |Vuj nj j_] as in the proof of this lemma, this yields: 


|^x,yVinJ/ | d~ | 

< K[C e u,CZ e9 ](^\Vx, y Pin\+l 2 \Vl, y Pi n \+'Y 3 \V 3 Pin\ 


+ I U L,// l + IVrc.yU*,,^ 


N 2 x,y U as,//\ 


We obtain then: 



| ^x^y^inj! | T |Vt7n,-L|~ 


< c 


llu 


* 

as,// ’ 


H\B{0,L))f 



( 86 ) 


and we bound the last integrals on the right-hand side by computing pi n with respect to qi n and 
applying Propositions [lDl [Tl] and [12] as in the proof of Lemma [20] Note that g™ has only one 
component because div x . y u* s 7/ = 0. □ 

We end this section by considering the case where the aperture admits a cylindrical invariance: 
7 = 7 (-\/a: 2 + y 2 ). In this case we introduce (u 0 ,po) and ( u iiPi) the solutions to the Stokes system 
associated with boundary conditions: 


u o = u±(0)e z , 

u{ = u* (0) + (xd x u* ± (0) + yd v u* ± (0))e z , 
We note that straightforward computations entail: 



■ nder = 0, 


u*er(3Si), Vie {o,i}. 


So, we have indeed existence and uniqueness of the pairs (ui,pi)j = o.i- We also introduce vo = 
v[ug], vi = v[u{] the respective approximations of u 0 and Ui constructed applying the steps 
depicted in Section 13.11 We note that, due to the linearity of the Stokes problem and of our 
construction, we have: 

Uas = U 0 + Ui , V as = V 0 + Vi . 

First, remarking that u}(0) • e z vanishes (so that there is no logarithmic terms yielding from 
the application of Proposition 1101) , we reproduce the computations in the proof of the previous 
proposition and obtain at first: 


31 











Proposition 22. If h £ (0,1], there exists a constant K depending on dll, S , Cf / 9 , C cvx , C e u and 
L such that there holds: 


it 


|Vm 


Vvi | 2 < K\\u*-H 3 (B(0, L))\\ 2 . 


We complete the study of u as — v as by computing the asymptotics of the most singular term: 
u 0 — v 0 . We prove: 

Proposition 23. If h £ (0,1] and 7 is radial, there exists a constant K depending on dfl, S, C^ e9 , 
C C vx , Ceil and L such that there holds: 

J' |Vu 0 - Vv 0 | 2 < K\\u*-H\B{ 0 , L ))|| 2 . 

Proof. As the mapping u* 1 — y Uq is continuous H 2 (B(0,L)) —> H 3 (B(0,L)) D Hz(d£l), we treat 
only the case Ug = e z and we drop superfluous index 0. 

I 11 this cylindrical case, a particular feature of Vo is that, in the problem solved by qi„, the weight 
7 is invariant by rotation around the origin and the source term is a constant. Consequently, qi n 
is a radial function and an explicit formula is available (as in the proof of Proposition [15]). Up to 
assume that xl is also radial, we get that v iT1// is directed along the radial unit vector e r . More 
precisely, we have: 

Vin,//(x,y,z)=v r (r,z)e r , a.e. in Q L . (87) 

with an explicit formula in the aperture: 

3 r 

v r (r, z ) = ~^Tj{z - (h + 7 i 0 )))(- - 76 W), a.e. in Q L/2 . 

Following the proof of Proposition [21] in the previous section, to bound the distance between 
vq and uq, we compute integrals 


hn [w] 



V Pin') ■ W. 


Again, as in the previous section, explicit computations and integrations by parts yield that: 



*x,y v in,// 


w„ 



^ L ’ ^x,yUJj_ 



d z Vin t ±d z W±. 


Here, we introduce that 


d z W± = —diVa^W// d z V in ,X = —div x ,yVinj/ . 


( 88 ) 


We plug these identities in the above integrals and integrate by parts. Because of the radial form 
of v mj// (see (1871) '). we have 

^x,ydrV x ,yVin,// = ^x,y^in,// * 


This entails: 


For the last integral, we bound \^7x,yVin,i.\ similarly as in the proof of Lemma [20l Introducing the 
bounds on 76 , 7 t and 7 , this entails: 


x,yUin,l. \ 



I) 
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so that, introducing now that Pi n = XLQin and applying Proposition [IT] 



x,y 


Wl 2 <K{C r 2 e9 ,C eU ,C cvx } 



Applying Propositions flOl and H2l to q in entails: 



Vin,±\ <K[C r 3 e9 ,C ell ,C c , 


Then, we truncate w 77 with Xl /2 that vanish on the lateral boundaries of <9<7 l/2 (and is equal 
to 1 on Gl /4 ) and we obtain that: 


>Ql 


^x,yVin,// * W// 


< 


'Ql /2 


|A x> yVi n>// \- |Xl/2W /7 | +C L \\v inJ/ \H 1 (VL\g L/i )\\\\™ // ]H 1 (yL)\\. 


With similar arguments as in the proof of Lemma [lH] (see (15U1) ), we bound the last term on the 
right-hand side: 

\\v inJ/ ;H 1 (n\g L/4 )\\ < K[C ell ,C r 2 e9 ,5}. 

Finally, we apply the formula for Vin,// which implies in particular that: 


^x,yVin,// — 


rdrr + 3d r (^) 


Hence introducing that w vanishes on the upper and lower boundaries of g ^/ 2 , we bound with 
Cauchy Schwartz inequalities and a Hardy inequality in the z direction: 


!■ p2 tt pL ph+-yt (r) 

/ |A Xi 3 / Vj „ i77 1 • |xl / 2 w 7/ | < / / / (z - 7 fc) IA 

JQl do do 


*x,y v in,// | 


|X-L/2W 77 | 

Z 


drdzd(9 


C 


' Ql 


7 


rdrr (^) + 3d r (^) 


|Vw;L 2 (J-)||. 


We emphasize that the constant C is universal and in particular independent of h. With the 
explicit formula for v r we get: 


7 


As 


rd rr (jL^j + 3 dr (y)] | < K[C r 2 e9 , Ceil} 

L ' 2 r 3 dr 


' Ql /2 


< c 


7(r) 




we get that, for a constant I\ depending on C^ es and C e u, C cvx , there holds: 

\Ii»M < iX||u*;iJ 2 (H(0,L))||||Vw;L 2 (J')||. 


This ends the proof. 


□ 


4 Proof of Theorem [5] 

In this last section, we exhibit a particular case where the above informations yield a sharp 
asymptotic expansion of the quantity: 

X |vu|2 ' 
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Throughout this last section, we assume that fl = K 3 \ B* with B* a sphere of radius R and that 
B is a sphere of radius S. We recall that we have then: 

7 t(x,y) = X + S U +0{(x 2 +y 2 ) 2 ), 

lb{x,y) = ~ X + 0((x 2 + V 2 ) 2 ), 

7(*,y) = ^-^y^ + o^ + y 8 ) 3 ), 


2 i?! 


8i?l 


where R\ and R 3 satisfy (12211 . We fix also a smooth boundary data u* and denote by u the exact 
solution to the Stokes problem with boundary condition u* and v[u*] the approximation that is 
constructed in Section 1371 

We introduce the notations of the previous section: indices 0,1, as, R distinguish the compo¬ 
nents of u* and u and v[u*]. Hence, we have 

v[u*] = v 0 + Vi + v R , 

and a similar decomposition for u. We decompose 

Il u - v[u*]; V\\ < ||u 0 - v 0 ; V\\ + ||ui - vi; V\\ + ||u fl ; V\\ + ||v fl ; V\\ 

Applying Proposition [T9l and Proposition [22l Proposition [23] in the cylindrical case we deduce: 

j f |V(u — v[u*])| 2 = 0(1) {||u*;^ 3 (H(0,L))|| 2 + ||u*; ^(9H)|| 2 } , 

where we keep the convention that landau notations 0(1) stand for quantities depending on 
(h, R, S ) which remains bounded by a constant depending on R , S only for h £ (0,1]. 

The weak formulation of the Stokes problem (remarking that u — v[u*] vanishes on dT) yields 
also that: 

J Vu : V(u — v[u*]|) = 0 , 

Hence, we have: 

j |Vu| 2 = J |Vv[u*]| 2 + 0(1) {||u*; H 3 (B(0, T))|| 2 + ||u*; ^(5H)|| 2 } . 

To compute the first integral on the right-hand side of this last equality, we split: 

[ |Vv[u*]| 2 = ^ Ei + 2 (E 01 + E or + E 1R ) , 

where for (i.j) £ {0, l,i?} we define: 

Ei= [ |Vv,| 2 


Eij = I Vv, ; : Vv,. 


IT 


It 


We complete the proof by studying the asymptotics of all these integrals. The first-order terms 
will yield by computing Eq and E\. 

4.1 Study of positive terms 

We recall that, by applying Proposition [El we get at first that 


I Vvr| 2 = 0(1) 1 1| u* ; H 3 (B(0, T))|| 2 + ||u*; (9S1)|| 2 | . 


(89) 
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Concerning the other terms, we remark that, by construction Vj = v.;.^ + v * >ex t so that: 


|VvJ 2 = / |Vvi,i„r + 2 / Vv,j„ : Vviext + / IVv,. 


/.f 


/F 


where, reproducing the computations in the proof of Lemma fl8l (see (ISHll h we obtain that, for 
i = 0,1: 

f |Vv liCTt | 2 + f |Vv M „| 2 = O(l)||u* ;j ff 3 (B(0,L))|| 2 . 

Jo* Jn* 

Hence, we get that, for i = 0,1: 

Ei = E i}in + O(l)\\u*-,H 3 (B(0,L))\\ 2 , Ei, in = |Vv^„| 2 . 

4.1.1 Asymptotics of £d,m 

Let first consider Vi. We drop index 1 in the sequel and we recall that vi j n , denoted here by v,; n , 
is constructed as follows: 

Vin,//{x,y,z) = ^(z-{h + 'Yt))(z-'r b )V Xt yPin- Xl u */ (0), 


= 2 diV:r ^ 
h+7t 


rh+ 7t 

/ (s- (/i + 7t))( s -7fo) v ^,2/P™ds 

J Z 

s- {h + 7 1 ) 


div 




7 


Xlu* (0)ds 


for (x,y,z) € where: 

Pin{x,y) = Xl(x, y)qin(x, y) , V (x, 2 /,z) e <7z, ■ 
with xl a suitable truncation function and qi n the unique solution to 


-^div a; , y ( 7 3 V a; ,j / gi n ) = (xd x u* ± (0) + yd v u* ± (0))e z - div Xjy 

Qin — 0 ; 

From Proposition II11 we obtain first that: 




on 5(0, L ), 
on dB( 0, L). 


Ei= ( |Vv m | 2 + O(l)||u*;f7 3 (l?(0,L))|| 5 

7Cr.,, 


>Gl, 

hence we may replace pi n by qi n in computations from now on. Then, it comes from the proof of 
Proposition [5l] (see (1851) 1 that: 


We obtain: 


[ |V„, y v mi// | 2 + |Vu m ,jJ 2 = 0(1)||u*; H 3 (B(0, L ))|| 2 . 

JGl/2 

E 1 , in = f \d z v lnW/ \ 2 + O(l)\\u*-H 3 (B(0,L))\\ 2 . 
JGl/2 


(90) 


Explicit computations yield that, on Ql/ 2 i there holds: 


dzVinj, = 


(z — (h + 7t)) + (z — 7&) ^ u //(0) 

q V x,y(lin 
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Consequently, we have that: Ei tin = E\ in + Ef in where (note that the cross-term vanishes by 
integrating w.r.t. ^-variable at first) 


E l,in = 


1 


/ l(~-(^ + 7t)) + (~-7b)| 2 |V*, y (M 2 > E 1 ,in= f 

JQt./ 1 ?. JQi 


| u *( 0)| 2 


*Gl /2 JQl /2 'y 

We end up the proof by computing the asymptotics of E\ in and Ef in . 

Concerning E\ in , we expand, for sufficiently small rg : 

r L / 2 r( \ r 

El,in = l ^ 


- 2 «(f ^ + Od,) 

= 2tt j R 1 |u*(0)| 2 | ln(fc)| + 0(l)||u*; H 3 (B(0, L))\\ 2 . 


As for E\ ln . we go back to the computations of Section [2] Indeed, integrating at first with 
respect to z, we get: 

1 / 7 3 |Vg in | 2 , 

Ib(0,L/2) 


E l,in 12 


where we apply Proposition[l6]to qi n to compute the asymptotics of this last quantity. This yields: 
Kin = ^Ri\Hh)\\RiV x , v ul(0)+ 2 ( ( ^~^ ) ) u // (0)l 2 + 0(l)l|u*;g 3 (i?(0,T))|| 2 . 


Finally, we obtain: 


24tt 


E 1m = 2irRi |u*(0)| 2 + —R^R^^yuK 0) + 


(R-S) 


-//v-n ' 5 1 2(i? + 5)“ // 


u //( 0)| 2 ) | hi(/i)| 


O(l)\\u*-H\B(0,L))r. 


4.1.2 Asymptotics of E(),in 

We focus now on E 0 , in and drop index 0 for simplicity. Let first recall that v in is constructed as: 
Vin,//(x,y,z) = 7}(z~ (h + y t ))(z-'y b )V x , y Pin, 

Vin,±(x,y,z) = idi v x ,y 
for ( x , y,z) S Ql, where: 

Pin{x , y ) = xl{x, y)q in {x , y), V ( x , y, z) &Q L . 
with Xi a suitable truncation function and qi n the unique solution to 

-^ div (7 3 Vgm) = M l(°) on .8(0, L), 

qin = 0, ondB(0,L). 

We recall further that, in this radial case, we may compute qi n explicitly: 

r L 6 s 

q in (r) = u* ± (0) -^rrds, Vre(0,i). 

Jr 7 


fh+ 7t 


(s - (h + 7t))(s - 7 b)^x,yPin ds 
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so that 


Vr G (0, L). 


d r qin(r) = ~u* ± (0) 

This entails that v^ n = v r e r + v z e z with: 


6 r 

7 3 (r) ’ 


v r (x,y,z) = 


d r C[in 


{(z - (/i + 7t(r)))(2 ~7b(r))) , 


= (( z ~(h + 7 t(r)))(z - 7b(r))) 


V (x,y,z) S t/ L / 2 ■ 


We assume from now on that m(]_( 0) = 1 for simplicity. We note that we have 

K(0)|<O(l)||u^tf 2 (B(0,A))||, 

so that all 0(1) terms in the following computations will turn into 0(l)||u*; H 2 (B( 0, L ))|| 2 in the 
final result. 

As in the computations for ifiyn, from Proposition II 11 we obtain first that: 


-E'O ,in — [ |Vvi„| 2 + 0(1), 

''Gl/2 

and we replace pi n by qi n . Also, we already computed in the proof of Proposition E51 that 

J\Vx,yVin,±\ 2 = 0 ( 1 ). (91) 

Consequently, we have: 

+ \d z V r \ 2 + IdzVzf) + 0(1), 


-S'O ,in — / ( |^r^r| 

''Gl/2 \ 

where, due to the incompressibility condition satisfied by v* 

2 


*Gl/2 


\d z v z 


' Ol/2 


d r \rv r 


Gl/2 


l ('* 

J yjL /2 \ 


\d r 


0 ( 1 ) 


/ / 


z=ft+ 7 t(L) 


—L J z=7b(L) 


as v r remains bounded independently of h away from the origin. We get thus: 

Eo,in = A + 2/ r + 0(1), 

with: 


A = 


'Sh, 


15, 


A = 


/a 




/a 


|<A 


Computation of I~ Replacing the integrand in A with its values yields: 

^ r fh+Jt(r) 


A = 


4 Jb( 0,L/2) J lb (r) 

1 

/B(0,L/2) 


|5 r gin(r )| 2 [(z - (/i + 7 t (r))) + (2 - 7 &(r))]“ dzdzd?/ 

|Vg in | 2 7 3 dxdy f ( 2 s - l) 2 ds 
Jo 


1 

12 


Ib(0,L/2) 


|Vg m | 2 7 3 da;d2/ 
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At this point, we apply Proposition [T5] to qi n yielding: 


/B(0,L/2) 


|Vg in | 7' dxdy 


= 72tt 


e> 2 p4 


0 ( 1 ). 


This entails finally: 


4 = - 18tt^|| ln 4)l + o(i) ■ 


Computation of I r We proceed with the computation of I r = I* + I* where: 

2 


Ir ■■= / \drVr \ 2 , 

JGl/2 

We first compute 7 2 by replacing v r with its values: 


It := 


[ 

v r 

’Gl/2 

r 


e = r p< r > [(^-(>.+ ^m,-^rw rdrdz 

Jo 4 


= 187T 


' o bi( r ) 

^/ 2 rdr f 1 , 


7(r) e 


[ f [(s-l)s}- 

Jo 74) Jo 


o 7(4 JO 
37t /' l / 2 rdr 

~5~ 


ds 


/o 74) ' 

We already computed (see the computation of Ef in ) that: 

ri/2 rdr 

A) 74) 

so that we obtain finally: 


= i?i|ln(/i)|+0(1). 


‘X'TT 

P r =^ Rl \\n{h)\+0{l). 


Second, we expand d r v r . We have: 


1 


d r v r = - ( d rr q in (r)(z - 4 + 714)))4 - 7b 4)) + 4<7m4)4[4 - 4 + 7t4)))4 “ 7b 4))]) 


Consequently: 


7r 4 L 


\d rr qin(r){z ~(h + 7t4)))4 - 7&4))| 2 

Gl/2 

+ f \d r qin(r)d r [(z - 4 + 7*4)))4 ~ 7&4))]l 

^Ql/2 


+2 / d rr q in d r q in d r [(z - 4 + 7t4)))4 “ 7*>4))]4 - 4 + 7*4))) 4 - 764)) 
JGl/2 

= Ir,a + Ir,b Ir,c 

After tedious but straightforward computations, we get: 

4,a = 157 t4i| ln(/i)| + 0(1), 

4,6 = ( 247^-^^11441 + 0 ( 1 ), 

A,c = —307ri?i| ln(/i)| + 0(1). 
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and 


It = — |ln(ft)|+0(l). 


Combining computations of I r and I z , we obtain: 



4.2 Asymptotics of cross terms 

We proceed with the computation of the asymptotics of cross terms , i ^ j € {0,1,-R} 2 . We 
recall that with similar arguments as previously, there holds: 


Eij = E^ in + 0(1) |||u*; H 3 (B(0, L))|| 2 + ||u*;if*(Sfi)|| 2 } , 


with obvious notations. 

4.2.1 Asymptotics of Boi.m 

We first treat the term Eqi. For this term, we assume without restricting the generality that 
l?Vuj_(0) + u,* (0) is parallel to ei. As a consequence, we obtain that q\j n reads ip(r) cos(0) in 
polar coordinates so that, as a function of (x, y) it satisfies: 

Qi,in(~ x > -y) = -Qi,in( x > y ) v ( x > v) e l). 

Hence, introducing the rotation matrix, 



there holds: 


SiVi,in{Si{x,y,z)) = -vi ,i n (x,y,z ), \/(x,y,z) e Q L 


On the opposite, we have that qo t in satisfies qo,i n {~x, — y) = qo,i n (x,y) so that we have the 
symmetries: 

Siv 0}in (Si(x, y, z)) = vo,in(x, y,z ), V (x, y, z) £Q L . 

Going to the derivatives, this entails that: 



4.2.2 Asymptotics of EoR,i n 
By definition, we have : 



We introduce: 


PR,in — PR,in 


/ PR,in{x, y)vo,±,in(x, y, Jb(x, y))dxdy, 
Jb(o,l) 


and, applying that VRj n is divergence free, we transform: 


/ - PRjnh) ■ Vv 0 ,in 

JQl 




V Rjn, _L ' V x,yVo,in,A. 
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In this identity, we integrate by parts the integrals on the first line of the right-hand side yielding: 

I x,yVR^in,// PR,ivX 2 ) • ^x,y^O,in,// “ 1 “ I ^z^R,in,// ’ ^z^Q,in,// 

JQl JGl 

= ^ (^x,yVR : in,// V x,ypR,in ) * Vo, in,// H" &zzVR,in,// * Vo,m,//^ 

= {^j i^x,y^R,in,// ^ x,yPR,in ) ‘ Vo, in,// “I” ^zz^ R,in,// * Vo,m,//^ 

where we used that 

• v o,m,// vanishes on 8Ql 

• ^zzV R y in,// V ’ x,yPR,in = 0 in 

We also note that pu : i n does not depend on z and apply boundary conditions for vg,±,i n yielding: 


)Ql 


iWm^A0,m,_L — 


I PR,in^O,in,l. (x,y,ib(x,y)) 

Jb(o,l) 


= 0. 


because of our choice for pR t i n . Finally, we have: 

Ao R,in = I ^x,y^R,in,// ' Vq, m,// A / ^z'^R,in ,.L ' ^z^0,m,_L 

JQl JQl 

A / ^x,yVR.in, A. ' ^x,y^0,m,J_ 

JQl 

From (l89l) and (l9ll) . we have that: 

J \Vx, v vo,in,x\ 2 A J' |Vv m ,*| 2 = 0(l)(||u*; H 3 (B(0, L)) II 2 ). 
Consequently, it remains: 

-^0 R,in = / ^x,y^R,in,// ’ Vo,in,// A I dzVR,in,A- ' ^z^0,m,-L 

•/ez, -tei, 

+0(l)||u* ; /r 3 (B(0, i))|| 2 

Applying that v/jy ri and vo,j n are incompressible, we get: 

j d z Vfi^ n ^d z VQ^i n ^j_ = / div^^Voyn,// diVx^Vfl,^,// 

■/ch ./.f 

As Vo ,//,in vanishes on 5C/£, we integrate this identity by parts: leading to: 

I ^x,y~^R,in,// ’ VOjin,// A / ^z^i?,m,J_ ‘ ^z^0,m,J_ 

Jg L -Igl 

~ I (^x,yVR y in,// A Vx,ydiVx,y ^R,in : // ) 'Vo, 

JQl 

so that forgetting remainder terms for conciseness: 

\EoR,in\ A / l^ 7 x,y V -R,m,// ll v 0,in,// |- 


™,// 
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At this point, we recall the explicit form of VR t // t i n : 

=(z-(h + 7 t ))(z - 7b)VPR,in + Xr."*n.,/ ■ 

Consequently, applying (1531) there holds: 


\^ x ,y^R,in,// I < K[R, S\ l\^PR,in\ + 7 2 \^'"PR,in\ + 7 I ^ PR, 


„ |Vub .. I |u* B .. I 

+ |V 2 uJj >7/ I + --M 


7 2 


7 


We recall we have also: 


1^0,in,// | < K(0)|-. 


Finally, applying that r < K[R,S]j3 , we get: 


EqR, in < A'[-R,S']|ul(0)| f (7"|Vpj?,i„| +7 2 |V 2 p fl ,in| +7"|V 3 p fl ,in|) rdr 


V 2 u* R j/ \ + \Vu* R j/ \ 


l R,ll 


rdr. 


Hence: 


EIr M < if[i?,5]| U l(0)|| I J (7 3 |Vp«,in| 2 +7 4 |VVi„| 2 +7 5 |VW| 2 ) rdr 


+ ll u ij //; H 2 (B(0, L))|| >, 


so that, introducing qR^ n = q^\ n + hq^\ n and applying Propositions l9l and 1T21 we obtain as for 

(ESI): 

E 1 0RAn =O(l)\\u*;H 3 (B(0,L))\\ 2 . 

4.2.3 Asymptotics of EiR^ n 

The computations of E\R^ n follows the line of the preceding section. We first remark that 
and ((90)) imply: 

J \V x , y vi,in\ 2 + |V^,i„,x| 2 + j r |Vv*, m | 2 = 0(l)||u*; H 3 {B(0, L))\\ 2 


so that: 


EiR. in = J d z v ltinW/ d z v RMt/ , + 0{1 )\\u*;H 3 (B(0,L))\\ 2 . 


Explicit formulas yield that, for i = l,R, we have: 


d z Vi t in,// = V pi tin {{z - (h + 7 1 )) + (z- 7 b)) + 
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Hence, integrating at first w.r.t. z and deleting vanishing terms, we obtain: 


dzVl,in,l/dzVR,i n ,f 


>T 


< K[R,S] 


IB(0,L) 


\^PR,in\\^Pi,inh 3 + 


/ B(0,L) 


< K[R,S](f |V Pl , m | 2 7 5 +\S7 PR . in \ 2 r- 

VB(o,L) 


l u k//l 2 , l u v/l 2 


Jb ( o , l ) 72 72 

As we already computed several times, the vanishing properties of u*# in the origin imply that 


|u* D „ | 2 lut 12 


1 R,// I 


I - 3 -b -—< 0(1)||u; H 2 (B(0, L))\\' 

Ib(o,l) 7 2 7 2 


and, applying Proposition [ill and fTOl we get: 

/ |V 9 l ,i n | 2 7 5 + |V te , m |V < K[R, S]\\u*;H 3 (B(0, L))f 

Jb(o.l) 


Hence, we have finally: Ein tin = 0(1)||u; H 3 (B(0, L))\\ 2 . 
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A Proof of Lemma fT7 

This appendix is devoted to the proof of 

Lemma 24. Given R > 0, there exists a q £ C oo ((0, oo)) 0 C([0, oo)) solution to 


d ss q 


3 s 


, s W + h) 


d s q-^ = - 
s z 


12s 

s £ (0,oo), 

(92) 

(i + m) 3 ’ 

lim q(s) = 0 . 


(93) 


Furthermore we have the asymptotic description: 


q(s) = 


48 R 3 
5 s 3 


1 144 R 3 1 

+ 0(~ j) d s q(s) = -^+0(- K ). 


(94) 


Proof. From now on, we let R > 0 as in the statement of our lemma. As the proof seems standard, 
we only sketch the main steps. 

First, we fix L > 0 and solve: 


12 dssQ 


1 


3s 


, s R ( 1+ 1r) 


dsQ~^ = 


(1 


2 R> 


s e (o ,L) 


9(0) =0, q{L) = 0. 


(95) 

(96) 
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To this end, we introduce 71 (x,y) := (1 + ( x 2 + y 2 )/(2R)) and introduce the bilinear form: 

((®,<?))i = To I |7i| 3 Vra7 : Vip 
Jb(o,l) 

on Hq(B(0, L)). As 71 is smooth and does not vanish on B(0,L) we obtain, by applying the 
Stampacchia theorem, existence of a unique vjl G Hq(B(0, L)) solution to 

((®L,P)) i=[ xip(x,y)dxdy, V <p € H%(B(0, L)). 

Jb(o,l) 

Due to the invariance by rotation of 71 and the computation domain B{ 0, L) for this weak formu¬ 
lation, we have that the unique solution reads in polar coordinates: 


uj l {s, 0) = q L {s) cos(0) V (s, 0) G (0, L) x (— 7 r, n) , 

with ql G C°°((0,L)) nC([0,L]) solution to (I§5l) - (1M1) . 

Second, we remark that we have a maximum principle for (1951) and that 

(2RY 


(1 + H) 3 ^ 


Vs > 0, 


this yields in particular that 

0 < q L (s) < 12(2R)i , VsG(0,L), VL>0. 

Furthermore, setting q(s) = 48i? 3 /(5s 3 ) there holds: 


1 

12 


&ssQ ' 


1 


3s 


, s w + h) 


where 


|e(s)| < 


Co 


d s q - 


Vs G (1, 00 ), 


(1 + iS) 3 


e(s) 


Setting also q(s) = 1/s 4 , there exists a constant c > 0 such that: 


-4 


1 


3s 


*(1 + h) 


d s q~ 4 I > 4- Vs >> 1. 


Finally, by application of the maximum principle, we obtain that there exists a constant K inde¬ 
pendent of L for which: 

18 A 3 K 

qL{s) = - 7 ^ 3 - + qrem(s) with \q r em{s)\ < ■ 

As the previous estimates are independent of L we can pass to the limit in L —> 00 and obtain 
in the limit a solution to (l92l) - (l93l) with the expected zero-order asymptotic expansion. We do 
not detail further this passage to the limit. As 71 remains strictly positive on (0,oo) we may 
apply classical cllipticity results to yield that the solution-limit is indeed smooth on ( 0 , 00) and 
continuous in 0 . 

Finally, we remark that q' = d s q satisfies: 

+ «V(2B))V(.)) = gMa + ^/Cfl)) 3 _ 12 , 

s 

We obtain the expected asymptotic behavior of q' by integrating this equation between si = 1 
and S 2 = s > 1 . □ 
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